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We investigate the properties of electronic states and optical transitions in hexagonal GaAs quan- 
tum dots within Alo.3Gao.7As nanowires. Such dots are particularly interesting due to their high 
degree of symmetry. A streamlined postsymmetrization technique based on class operators (PTCO) 
is developed which enables one to benefit in one run from the insight brought by the Maximal 
symmetrization and reduction of fields (MSRF) approach reported by Dalessi et al. , after having 
solved the Schrodinger equation. Definite advantages of the PTCO are that it does not require 
having to modify any existing code for the calculation of the electronic structure, and that it allows 
to numerically test for elevated symmetries. We show in the frame of a 4-band k ■ p model that 
despite the fact that the D6h symmetry of the nanostructure is broken at the microscopic level by 
the underlying Zinc Blende crystal structure, the effect is quite small. Most of the particularities 
of the electronic states and their optical emission can be understood by symmetry elevation to D6h 
and the presence of approximate azimuthal and radial quantum numbers. 

PACS numbers: 42.81. Qb, 81.07.Gf, 42.55.Px, 



I. INTRODUCTION 

Semiconductor nanowires have emerged as promising 
building blocks for realization of various nanoscale op- 
toelectronic devices 2 . The nanowire technologies en- 
able heterostructures with a high level of flexibility in 
terms of geometry and material composition 2 . In par- 
ticular, growth of site controlled quantum dots (QDs) 
within nanowires show significant advantages compared 
to more conventional self-assembled QDs (Stranski- 
Krastanov) 3,4 , and bright single photon emitters have 
been demonstrated in a range of materials using nanowire 
QDs 37 . 

QDs may be used to facilitate quantum information 
and cryptography technologies, e.g. by emission of en- 
tangled photon pairs. In this respect nanowire QDs are 
especially suitable as they are highly symmetric, with a 
hexagonal cross section. High symmetry is a requisite for 
entanglement, needed in order to limit the fine structure 
splitting of the excitonic states which is a consequence 
of symmetry breaking. QDs within nanowires grown in 
the high symmetry direction [111] are therefore suggested 
as ideal sources for emission of entangled photon pairs 8 . 
Cascaded emission spectra, possibly enabling such pairs 
have also been measured experimentally 7 . Generation of 
entangled photons is also possible in self-assembled QDs, 
especially InGaAs/GaAs QDs grown on [111] substrates 
are proposed ideal for such generation, due to the three- 
fold rotation symmetry of this surface 9 . QDs may also 
enable quantum computation by using the electron spin 
as a quantum bit 10 , QDs grown in inverted pyramids on 
[111] GaAs substrates 11 are also promising in this re- 
spect. 



With increasing control of heterostructure shape and 
material composition, comes the possibility to grow 
structures according to optimized design. Theoretical 
models e.g. providing insight into the nature of electronic 
states are essential in such optimization. In this work 
we calculate the electronic states and optical transitions 
in hexagonal GaAs/AlGaAs nanowire QDs, and analyze 
the results taking advantage of their high symmetry. A 
detailed procedure for full-depth symmetry analysis and 
reduction of computational domain has been presented 
previously by Dalessi et al. 1 . Here we will pursue this 
effort by developing a new procedure alleviating recod- 
ing: a postsymmetrization technique using class opera- 
tors (PTCO). The PTCO is very general and applies in- 
dependently of the method employed for the calculation 
of the electronic structure (k • p, tight-binding e.t.c). It 
also provides a systematic and flexible procedure to test 
possible elevated symmetries. 

In ?? relevant theoretical studies of similar semicon- 
ductor heterostructures are considered, and we present 
the general features of the k • p model used in this pa- 
per. The explicit QD under consideration is described 
in Eq. (II. 3), and the symmetries of our model structure 
are identified, providing the premises to optimally choose 
the basis of our k • p Hamiltonian. Having fully specified 
the numerical model, we describe the pertaining sym- 
metry implications on the QD eigenstates in Eq. (III. 8). 
The PTCO is then presented in Eq. (IV.5). Eq. (V.13) 
and footnote b contain the numerical calculations, ana- 
lyzed using the PTCO. We show in Eq. (V.13) that the 
analysis gives rise to a deeper understanding of the elec- 
tronic states, in particular of the level sequences. In foot- 
note b we investigate the fine structure of the spectrum 
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of squared momentum matrix elements. We prove that 
symmetry elevation to D@h and the existence of azimuthal 
and radial quantum numbers are necessary ingredients to 
explain many missing/weak transitions. 



II. QUANTUM DOT DESCRIPTION AND K ■ P 
MODEL 



Despite large activity within the experimental realiza- 
tion of nanowirc QDs, there has until now been much 
less attention towards numerical calculations of the elec- 
tronic states. Niquet et al. 12 did perform calculations of 
strained InAs/InP nanowire QDs, using a tight binding 
model. The optical transitions were given, and labeled 
using group symmetry. The QDs were however approx- 
imated as cylindrical in that work; cylindrically shaped 
QDs grown in the [111] direction of a wurtzite structure 
will inherit the C3V symmetry of the crystal. Of related 
interest is also the work of Zhang et al. 13 , considering 
excitons in nanowire QDs in the strained InGaN/GaN 
material system, using an effective mass approximation. 
Nanowires in the GaAs/AlGaAs material system consid- 
ered in the current paper were treated by Kishore et al. 
using the k • p model, but no calculations exists, to our 
knowledge, of nanowire QDs within this material system. 

The k • p theory, originally intended for the calcula- 
tion of band structure of crystalline solids, has also been 
widely used for calculation of band structures in het- 
erostructures including QDs. Large emphasis has been 
on the strained self assembled QDs 15-19 , also including 
calculations on GaAs/AlGaAs QDs 19 . 

QDs with the particular hexagonal shape considered 
here were studied numerically using the k • p model in 
Rcf. 15, no explicit usage was however made therein of 
the symmetry properties. Ref. 20 present a strictly qual- 
itative symmetry study of electronic states and optical 
transitions in hexagonal QDs. 

In the fe • p model one assumes a weak interaction, 
fe • p, between the crystal momentum and the electron 
momentum 21,22 . Using Bloch's theorem, the wavefunc- 
tion can be separated into a slowly varying envelope func- 
tion, and a rapidly varying part with the periodicity of 
the crystal lattice. Note that this approximation is only 
valid close to the zone center (k = 0), and it is not capa- 
ble of describing an interface between materials of differ- 
ent crystal structures, e.g. a transition from a Wurtzite 
to a Zinc Blende material. 

The symmetry preserving ability of the k • p model 
has been investigated numerically 23 , demonstrating that 
the true symmetry of any structure can be restored upon 
inclusion of enough bands and interface terms in the k • p 
model. Care should nonetheless be taken to distinguish 
between the symmetry of any simplified numerical model 
and the physical system. 

We shall use a simple k-p model to describe the states 
of a GaAs QD within an AlGaAs nanowire. The elec- 
trons of the conduction band is described using an effec- 



tive mass model, and the holes of the valence band are 
described using a 4-band Luttinger Hamiltonian. This 
rather simple model will be used to demonstrate the 
PTCO, which can also be used in more complex theo- 
retical frameworks. 

The conduction band describes electrons with spin 
j = 1/2, the electrons can however be written using a 
scalar Schrodinger equation with an effective mass ap- 
proximation, ignoring mixing with other bands 1 
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Here, V is the 3D differential operator V = ^u x + 
^u y + JjU z , m* (r) is the effective electron mass in units 
of the electron mass m , and V^ b (r) is the effective con- 
finement potential for electrons in the conduction band. 
The envelope function tp n of energy level n is given by 
the Schrodinger equation: 



(11.2) 



The spinorial nature of the conduction band states can 
be restored later; the exact procedure is given in 24 . 

Band mixing and spin cannot be ignored for the holes. 
The top six valence bands can be described as multiplct 
states with spin j = 3/2, and j = 1/2. 21 The latter 
multiplet is the split-off band, separated from the first 
multiplet with the amount A so . If coupling to this split- 
off band can be ignored, we are lead to a four valence 
band k • p envelope function model. We can use the 4x4 
envelope function Luttinger Hamiltonian for diamond, 
if we neglect in addition inversion symmetry breaking 
in GaAs/ AlGaAs. This can always be expressed in the 
form 21 ' 22 . 
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Here, p, q, r, s are second order polynomials of differen- 
tial operators. Their exact expressions depend on the 
Bloch basis which will be chosen later after considering 
the symmetry of the Hamiltonian. The polynomial co- 
efficients are given in terms of the Luttinger parameters 
7i(r), i = 1 . . . 3, and Vvb(r) is the confinement potential. 
The valence band spinors tp are found using Eqs. (II. 2) 
and (II. 3). The underline will be used throughout to 
distinguish spinors from scalar functions. 



III. MODEL STRUCTURE AND ITS 
SYMMETRY 

In this section we identify the symmetries of our model 
structure, introduce some concepts from group theory, 
and review its implications for the best Bloch basis. This 
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basis will be given from the Maximal symmctrization and 
reduction of fields (MSRF) technique, previously pre- 
sented by Dalcssi et al. 1 . 

The model structure considered in this paper is a QD 
grown as an axial insert of GaAs within an Alo.3Gao.7As 
nanowire. A radial shell of Alo.3Gao.7As is grown around 
the dot so that it is surrounded by Alo.3Gao.7As in all 
directions. The growth direction defining the nanowire 
axis is the crystal direction [111], and the cross-section 
is hexagonal. We will assume a Zinc Blende crystal 
structure. A schematic of the structure is shown in 
Fig. 1. Similar QD structures have been grown previ- 
ously by Kats et al. 5 , however, they did obtain mixed 
crystal phases containing both wurtzite and Zinc Blende. 
Pure Zinc Blende GaAs/AlGaAs axial heterostructures 
has however been realized 25 . 
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Crystal structure: Zincblende 

FIG. 1. Schematics of a GaAs quantum dot within an 
Alo.3Gao.7As nanowire 

The true QD symmetry is the common symmetry of 
the mesoscopic hcterostrocture and the microscopic crys- 
tal structure. The mesoscopic structure (Fig. 1) is sym- 
metric (i.e. invariant) under discrete ±27r/6 rotations 
and under mirror operations w.r.t the six vertical planes 
containing the rotation axis, as well as w.r.t. the hori- 
zontal plane orthogonal to the rotation axis. These sym- 
metry operations and their compositions form a point 
group of 24 elements called Dgh 26 - Mathematically, the 
mesoscopic symmetry restriction on the conduction band 
Hamiltonian is expressed by the invariance relation per- 
taining to the confinement potential and effective mass 
appearing in Eq. (II. 1): 
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Here, 1Z is a set of standard representation matrices op- 
erating in 3D space 1,27 , and TZ(g) represents change of 
coordinates indexed by the group element g. The sym- 
metry of our conduction band Hamiltonian (Eq. (II. 1)) 
is then directly given by the mesoscopic symmetry Dgh- 



For the valence band, there are similar invariance con- 
strictions due to the mesoscopic symmetry on V v b(r), 
and also on the spatially dependent Luttinger parame- 
ters 7i(r) (appearing in the p, q, r, s polynomials) 
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However, the valence band Luttinger Hamiltonian 
(Eq. (II. 3)) also carries a face centered cubic diamond 
Oh symmetry due to the underlying crystal, contained 
within the first term of Eq. (II. 3). When the orientation 
of the crystal axes w.r.t the mesostructure is as in Fig. 2, 
the common symmetry elements are those of the group 
C 3v , i.e. {e, C£, C3", a v \, a V 2, Uvz}- Here and C% are 
discrete 2-7r/3 rotations and a v i, i = 1 ... 3 are three verti- 
cal mirror operations shown in Fig. 2. Hence the valence 
band Hamiltonian (Eq. (II. 3)) has only C3V symmetry. 
In the following we will therefore use the C3V group as 
the main reference when introducing the necessary group 
theory concepts. 

Description of the group, including group multiplica- 
tion tables and group representations can be found in 
reference books like Altman 26 . For the single group C3V 
there are three irreducible representations (irreps), E, A\ 
and Ai- The E irrep is a 2D representation whilst the 
Ai irreps are ID. 

So far we have just considered the single group, de- 
scribing operations performed on the spatial coordinates. 
The spinorial nature of the conduction and valence band 
does however necessitate the use of a double group 
representation 28 . In the C3V double group there are two 
ID irreps 1 E 3 / 2 and 2 E 3 / 2 and one 2D irrep Ei/ 2 - 26 - 
With spin a change of the spatial coordinates i}(g)^ 3D ^ 
('&(g)( 3D,> is an abstract functional operation correspond- 
ing to the change of coordinates described by the matrix 
7Z(g)) should be accompanied by a corresponding change 
in spin space, i}(g)^\ Such composite operations arc 
thus defined by fl(g) = tf(g)( 3D ) ® i?(3) (i) , where <g> is 
the tensor product between both operator spaces, and 
g becomes an element of the double group. A natu- 
ral choice for a standard representation of rotations in 
spinorial space in terms of Eulcr angles, g = (a, f3, 7), is 
the set of 4x4 Wigner matrices, W(a, j3, 7) 27 , since the 
Luttinger Hamiltonian is usually expressed in a Bloch 
function basis transforming like angular momentum and 
indexed \j, m), j = 3/2, where m is the component along 
some quantization axis. Any matrix for improper rota- 
tion (i.e. mirror operations) can be obtained as a com- 
bination of proper rotations and spatial inversion, i, but 
since the diamond crystal structure is even under i, this 
can be ignored. 

To fully take advantage of the symmetry properties of 
the valence band Hamiltonian, a specific choice of the 
Bloch basis is needed 1 . For some low symmetry groups, 
an optimal quantization axis can be found 1 , enabling 
full depth symmetry analysis. However, as explained in 
Ref. 1, in general a rotation of the quantization axis is 
not sufficient, and a more general unitary transformation 
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of the basis has to be performed. The elements of the 
optimum basis are called 1 Heterostructure Symmetrized 
Block Functions (HSBF's), and are labeled by irreps of 
the double group. 

The HSBF's arc symmetrized superpositions of the 
usual Bloch functions. For our C3V Hamiltonian, the 
expressions for the HSBF's are given by Eqs. (B.5) 
and (B.6) of Eq. (A.6). The HSBF basis is written in 
terms of them as 



{| 1 £'3/2) 7 |-^l/2; l) 1 |-^l/2l 2 ) , \ 2 E 3 / 2 )}. 



(III.3) 



The ordering in Eq. (III. 3) is important to preserve the 
form of the time reversal symmetry operator 24 . 

Note that our choice of axes (Fig. 2) differs from that 
of Refs. 1,24, in which the z-axis in the equations for 
the HSBF's pertain to a quantization axis in the crystal 
direction [110], i.e. corresponding to the y-axis in our 
present choice of coordinates. We have rather chosen 
the z-axis based on the special axis of the nanostructure 
([111]), as the appropriate concept of light hole (LH) and 
heavy hole (HH) for our structure relates to m z 



±1 or 



±| respectively for a z-axis along [111]. 



In the HSBF basis, the p, q, r, s polynomials appearing 
in Eq. (II. 3) are given by 
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withe = These expressions have been obtained by 

using the bulk expressions, e.g. by taking the Luttinger 
Hamiltonian expressed in direction [111] in Ref. 22, then 
changing its basis to the HSBF given in Eq. (A.6), and 
finally by replacing kj ~i-§j f° r 3 — %■> Vi z - 

In the forthcoming analysis, it will be useful to describe 
the spinorial nature of the quantum states using the con- 
cept of light hole (LH) and heavy hole (HH) dominant 
character (because of bandmixing the HH and LH states 
are mixed) . Raw approximations for the effective masses 
parallel and perpendicular to the [111] axis can be read 
from Eqs. (III. 4a) and (III. 4b) (assuming r=s=0), and 
are given by 
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(III.6a) 
(III.6b) 



where the upper and lower sign applies to HH and LH 
respectively. For a normalized state |^>), the weight of 
the LH contribution is defined by 



W ^ 1J (^) ^ U(3/2,l/2|^>r + U<3/2,-l/2|-0}| 2 , (IH.7) 



and similarly for the HH, so w 
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1 . Using 



the basis change between the HSBF basis and the basis 
with quantization axis [111], it is easy to show that 



<' W=\(E 1/2 ,l\i>)\ 2 + \(E 1/2 ,2\ip)\ 



(LIU 



We thus note that LH/HH character w.r.t. [Ill] is con- 
veniently in one to one correspondence with Ei/ 2 / l E 3 / 2 
character in the C3V HSBF. 



IV. SYMMETRY IMPLICATIONS ON 
NANOWIRE QUANTUM DOT STATES 

The HSBF basis allows to maximally symmetrize the 
envelope functions appearing in the valence band spinors, 
and to give interpretation of the spinors as dominant 
product states. These possibilities, which are conse- 
quences of symmetry, arc explored in the forthcoming 
subsections. 



A. Symmetry of envelope functions 



FIG. 2. Coordinate axes and crystal directions for the hexag- 
onal QD. 



We first recall the transformation properties of the 
double group spinors, in order to study the properties of 
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envelope functions occurring in our nanowire QD. This is 
based on the MSRF technique 1 . The HSBF-transformed 
Wigner matrices are denoted by V B (g), as they are block 
matrices in this basis 1 . Indeed, one can write them as a 



E3/2 , of the irrep matrices 



direct sum, 1 E 3 / 2 © ^1/2 

of the C3V double group 1 . Therefore the spinor transfor- 
mation (passive point of view), under any double group 
operation g, can then be written as 



[±{r)\' = = Yy B (9)]jmng)- l r)] k - 

k 

(IV.l) 

A seminal consequence of symmetry in a system sub- 
jected to invariance relations like Eqs. (III. 1 ) and (III. 2) 
is that it is always possible 29 to label an eigenstate ip by 
an irrep, T, and a partner function index, fj. of the sym- 
metry group, such that its corresponding transformation 
law is also: 



%)^( r ) = E^ r (ff)]^€( r ) ( IV - 2 ) 

v=l 

where D T {g) is a set of representation matrices for T, 
and d T is the dimension of T. For D 1 (g) we shall use 
the matrices of Rcfs. 1,27, which correspond to a trans- 
posed multiplication table w.r.t. Ref. 26 since we use the 
passive point of view. 

It is clear that Eqs. (IV.l) and (IV. 2) strongly con- 
strain the envelope function shapes, and it was shown in 1 
that they can be uniquely decomposed into ultimately re- 
duced envelope function (UREF) components, transform- 
ing according to single group irreps. The classification of 
UREFs using single group labels could be used to re- 
duce the computational domain significantly as argued 
in Ref. 1. We will see in the present work that it is 
also very helpful in the postprocessing and analysis of 
the data. In our QD, the full valence band spinors can 
be labeled by C3V double group irreps, and expressed in 
terms of UREFs as in Eq. (IV.3) 
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The envelope functions have here been labeled in a sim- 
plified manner, keeping the double group labels of the 
global spinor and of the HSBF 24 implicit. The full labels 
should be restored for description of how the envelope 
functions of separate subequations in Eq. (IV.3) trans- 
form into each other. 

It should be pointed out that for QDs time reversal 
induces a unique mapping between Kramers degenerate 
spinors, in analogy to the k — > — k mapping of quantum 
wires 24 . The time reversal operator, K, is 24 K = FKq, 
where Kq, is a complex conjugate operator and Fjk = 
*k,(6-j)(-l) 3 ', ihk) = 1...4. The 2D E 1/2 irrep is self 
conjugated, and the ID irreps are mutually conjugated. 
Therefore, the pair of cigenstates belonging to a given 
eigenvalue is either the two partners of the self conjugated 
2D irrep (£^1/2) or the pair of mutually conjugated ID 



'3/2 



^3/2 



irreps ( 1 £ , 3 / 2 , 2 E 3/2 ). Accordingly, Kf_ 
and Kip El/2 = tp El/2 . It is easy to see from the form of K 
that the UREFs appearing in (IV. 3a) must then be equal 
to the corresponding functions in (IV. 3b). Similarly, for 
the E 1 / 2 irrep one obtains restrictions on (IV. 3c) and 
(rV.3d): 
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i = 1,2. (IV.4) 



is imaginary 



In the scalar approximation, C3V conduction band en- 
velope functions can be given single group labels directly 
and states are either non-degenerate, T = A i: i = 1 ... 2, 
or twice degenerate, T — E. Accounting for spin, time 
reversal symmetry implies Kramers degeneracy and the 
cigenstates of the conduction band are either two-fold or 
four-fold degenerate. This also holds when one considers 
the full Dgh symmetry of the conduction band. 

Having established in the present section the precise 
symmetry of every spinorial component using UREFs 
bearing single group irrep labels, we now propose to 
quantify the respective contributions of each irrep, as- 
suming that the full spinors are all normalized to unity, 



i.e. 







1. 



The set of UREFs appearing in a given set tf> (r), 
transforming into each other according to D T (g) (of 
Eq. (IV. 2)) are specified by ^r^I ■ Here, label the rel- 
evant HSBF block and (T a ,fi a ) every individual UREF. 
There are internal redundancies within the UREFs due 
to their own transformation properties. We can define 
the weight of a subset of symmetry T a within the block 
Tf, of the spinor r (assuming that Tb occurs only once 

in T) 30 : 
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(IV.5) 







r,r 



r,r a 



is independent of /!. In addition 

independent of fi a by symmetry 30 (generalized Wigner- 
Eckart theorem) ; therefore partner function indices \i and 
jj, a are not relevant. With the chosen normalization, one 
of course has the completeness ^2 Tb r d-p a w r, " ra = 1. 
The redundancy by dr a is because every UREF partner 
function contribute equal weight. 



One last comment regarding the DP GPS and its re- 
lationship with the HH/LH concept: as explained in 
Eq. (II. 3), in C3V the HH- weight equals the combined 
weight of the two £"3/2 spinor components. The valence 
band state of the previous example (typically our QD 
ground state) may therefore be classified as HH-like. As 
shown in Eq. (B.6) the two corresponding "spin states" 
are distinct conjugated DPGPS in C3V, and to the two 
partners of £3/2,9 DPGPS in Dgh. Similarly LH-like spin 
states will be associated to the partners of £1/2 DPGPS 
in C3 V and to those of £1/2, g DPGPS in D 6 h- 



B. Product states, dominant HSBF as DPGPS 

In a k • p model, the global symmetry of the eigen- 
states of the heterostructure can always be expanded into 
sums of products between symmetrized scalar envelope 
functions labeled by single group (UREFs), and a 
HSBF labeled by double group Often, when band 
mixing is not too large, there is one dominant term in 
the expansion. One can infer that the first factor de- 
scribes well the spatial probability distribution, and the 
second factor entails the symmetrized spin information 
encoded in the HSBF. In such a case we will denote the 
relevant HSBF the Discrete Point Group Pseudo Spin 
(DPGPS) 31 . This notation express the fact that one can 
add such spinor parts in the same way as ordinary spin, 
but using the Wigner-Eckart theorem for point groups. 
Expressing the valence and conduction band states as 
such product states is often useful e.g. when build- 
ing excitonic and biexcitonic complexes, and for sim- 
ple understanding of degeneracies and other properties 31 . 
As an example, consider an exciton with C3V symme- 
try. Assume that the symmetry of the main electron 
states contributing to this exciton can be written as a di- 
rect product r olcctron = A\ ® Ex/2- Similarly, assume 
that the symmetry of the hole states can be written 
as r holc = A 1 ® ( 1 £ 3 / 2 ® 2 £3/2)- The excitonic states 
can then be constructed by combining separately the 
spinorial parts and the envelope parts. The symme- 
try of the excitonic states is thus T cx = (Ai (g) A\) ® 
(E 1/2 ® (^3/2 2 E m )) = 2E. 

For the conduction band, writing the states as prod- 
uct states is straightforward. It is more complicated for 
the valence band, as band mixing may prevent a descrip- 
tion of the states in terms of individual DPGPS. In the 
HSBF basis it is however easy to test the weight of ev- 
ery UREF part appearing in Eq. (IV. 3) with the help of 
projection operators for scalar functions (on single group 
irreps) and calculating the weight (Eq. (IV.5)) of every 
individual component. We shall see that it often leads to 
interpretation of the nature of the eigenstates as prod- 
uct states of a dominant envelope function with a cor- 
responding DPGPS. For example an 1 E 3 / 2 state (IV. 3a) 
with the dominant weight in the first (^4i) component, 



V. POSTSYMMETRIZATION 

Postsymmetrization has been first demonstrated by 
Gallinet et al. using projection operators 30 . We will 
now present the PTCO, a novel systematic and stream- 
lined postsymmetrization method developed to disentan- 
gle and classify numerical eigenstates such that they can 
be given individual symmetry labels. The advantages of 
the PTCO are: 1) that it can be used in combination 
with possibly existing calculation code for the electronic 
structure, so that little additional numerical/theoretical 
work is necessary (it is general and not limited to k • p 
theories). 2) That it will allow to test and qualify pos- 
sible approximate elevated symmetries. We stress again 
that though these techniques are here exemplified start- 
ing from a simple 4 band k • p model, the PTCO can 
equally well be applied in the frame of more complex 
models. 

For postsymmetrization, one considers a subspace of 
the space spanned by the set of eigenstates of the Hamil- 
tonian. This subspace S of dimension d$ is obtained after 
selection of a few relevant eigenstates, {ip , j = 1 . . .ds}, 

which are numerical solutions of Eq. (II. 1) or Eq. (II. 3). 
It will be called the solution space. Usually we shall 
retain a finite set of neighboring eigenstates, from the 
ground state to a given maximum energy, e.g. g?s=20 in 
Eq. (V.13). We shall assume that the states ib , are not 

entirely as symmetric as they should be (due to numer- 
ical errors or due to the use of a non-symmetric grid). 
Postsymmetrization techniques will aim at finding the 
best symmetrized states that can be found in S, on the 
basis of the ideal symmetry properties that such linear 
combinations should satisfy. 

For any operator, F, there is a matrix representation 
in S, called 0(F), defined as 



ib Fib 
—3 



(V.I) 



corresponds to the product state T 



hole 



A X (g) 1 E; 



3/2- 



In the following subsections all operators, including the 
Hamiltonian H, the symmetry operators $(<?), the pro- 
jection operators Pj, and the class operators Ci, will be 
identified with their matrix representation in S. 

Note that in the general case numerical eigenstates do 
not exactly obey Eq. (IV. 2) and they can be mixed in 
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two ways. First, mixing can be due to fundamental de- 
generacy and the arbitrary choice of basis by the solver 
within every eigenspace. This effect can be strong but 
is not in violation of the symmetry of the (analytical) 
Schrodinger equation. Disentanglement is however nec- 
essary to enable the use of individual group labels corre- 
sponding to standard irreps D T (<?) and belonging partner 
function indices /i. Second, eigcnstates can be mixed due 
to numerical inaccuracies, in particular if the grid does 
not respect the symmetry Such mixing is always weak, 
when good convergence is achieved. If energy levels are 
closely spaced, symmetry breaking due to the grid may 
nevertheless couple states at nearby energies and with 
different symmetries. The goal of the PTCO, which is 
a global postsymmetrization method, is to disentangle 
both kinds of mixing simultaneously and it will be pre- 
sented after a short description of a more standard tech- 
nique based on projection operators. 

A. Postsymmetrization using projection operators 

A standard procedure to analyze and postsymmetrize 
the computed basis {tp.,j = l.-.d^} is to use the a set 

of projection operators, Pj which project any state onto 
the partner function fi of the irrep T. The projection 
operators are given by 29 

where \Q\ is the number of elements of the group Q. A 
matrix representation in S can be built using Eq. (V.l). 
The set of projection operators are orthogonal, and 
obeys 

E^ r = 1 - ( v - 3 ) 

They can uniquely decompose any state into 

The fact that the grid slightly breaks invariance will 
break the expected symmetry properties of the com- 
puted eigenstates P^Vv = $r,r'^ ( u, J u'V'Ji- To quantify the 
amount of symmetry breaking for the basis functions of 
each double group irrep, we define the weight of (r, fj.) in 
spinor -0 by 

w l (t) = \Kt\\ ■ (v- 5 ) 

For a perfectly symmetrized state i/j, we have w^(ip) £ 
{0,1}. 



The starting point for the standard postsymmetriza- 
tion procedure is to compute the actual set of numbers 
w^(ipj),j = 1 . ..ds for all (r, fi), which give information 
on the symmetry content of any state. This informa- 
tion is used to identify subspaces where states of different 
symmetry have been mixed. 

Then in a second step, symmetrized subspaces for all 
irreps can be built starting from a suitable set of rep- 
resentatives. Let k be in this set, and let w^(4>k) be 
larger than the corresponding weight of the other ir- 
reps and other partner functions. The projected state 

ip = Ppiftk/ \\Ppil>k\\ transforming like is then a 
good starting point for the symmetrized subspace, the 
non-projective part of state k is thrown away. For the 
ID irreps e.g. l E 3 / 2 this concludes the procedure, but for 
higher dimensional irreps like E 1 / 2 the remaining partner 
functions must be found to span fully every irrep sub- 
space. One can then use transfer operators T^ v ,fi ^ v, 

where T? v = j§j J2 ge g (° r (fiO^Cff)) (a generalization 
of Eq. (V.2)), and the property tpj^ = T^ip^ to generate 
the basis. 

Rather than the standard procedure, we propose below 
to use a global procedure, which do not throw away any 
information by selecting representatives. Here, all states 
will be taken into consideration in one step, rather than 
operating on every degenerate subspace sequentially The 
procedure is not cumbersome, and it is easy to use. It 
does allow to automatically generate well behaved sym- 
metrized eigenstates of the Hamiltonian. 

B. Postsymmetrization using class operators 

This procedure is based on the concept of commuting 
class operators, which are related to conjugation classes. 
The classes Ci of a group consist of elements that are 
conjugate to each other. Two group elements g\ and g 2 
of the group Q are conjugate if there exists an element g 
in Q such that g 2 = <?(<7i)<7 • 

Class operators are a set of operators {Ci} constructed 
by summing operators corresponding to each clement of 
a given class 32 

d = <%)• (V.6) 

They are closed under multiplication and commute both 
with each other and with all group operations. These 
three properties make them useful e.g. in relation to 
the sought decomposition of the solution space into or- 
thogonal symmetrized subspaces. The number of class 
operators equals the number of irreps of a group. This 
yields three classes for the C3 V group. One class consists 
of the identity, C e = {e}, one class of the three mirror 
operations, C a = {cr vi ,i = . . .}, and the two rotation op- 
erations constitute the final class, C3 = {C3 ,C^}. The 

— r 

sought symmetrized basis vectors ip of S are eigenfunc- 
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tions of the class operators. The class operator Ci has new basis). Identifying the /c-th symmetrized state as 

— r 

one of the sought symmetrized "eigenvectors" ip , we get 



eigenvalue A*, i.e. G'iip 1 = A^?/> r . This eigenvalue can 



be found directly from the character of the class, xf ; one 
has 32 



A 



(V.7) 



Here \Ci\ is the number of elements in the class Ci. 
Eq. (V.7) can be used to deduce the group symmetry 
labels for the eigenfunctions directly without having to 
employ the projection operator. 

The power of the class operators has been extensively 
described by Chen? 2 . He showed in particular that 
the set of class operators linked with a canonical sub- 
group chain can be used to form a complete set of com- 
muting operators (CSCO) for the group space. In our 
case, since we want a symmetrized basis in the solution 
space, we must add the Hamiltonian to have a CSCO, 
since the eigenvalue of the Hamiltonian will be sufficient 
to distinguish subspaces with identical symmetry. For 
our C3V case, only one canonical subgroup is relevant, 
C s = {e,o~ v \}. It will distinguish the partner functions 
belonging to the 2D irrep E 1 / 2 - To this end we add the 
operation i^er^i), which is a class operator for C s . The 
CSCO in solution space is then: 



Vi 



1> 



4^1^- The overbar stresses that in fact they 



are not anymore numerically perfect eigenvectors of H. 
Nevertheless, they are in a way more physical since part 
of their due symmetry, which was broken by numerics, 
has been restored by the procedure. The index I was in- 
troduced to obtain a unique labeling, and is here usually 
chosen so as to number states with identical symmetry 
(r, n) by order of increasing energy. The symmetrized 
states can be easily constructed with 



— jU,i 



ip r . ) ip .. 

— /■»>*/ —3 



(V.ll) 



Having obtained a symmetrized basis, we must now 
identify the three labels T,/j,,l for all states. The eigen- 
values of C can be discarded. It is much more infor- 
mative to recompute the average values of the original 
CSCO (Eq. (V.8)), which will provide the labels for the 
symmetrized states. First, 



— r / — r 



Hip 



■fi.l 



(V.12) 



CSCO = {H,C a ,C a ,4(a vl )}. 



(V.8) 



Here the trivial class operator C e has been omitted since 
it corresponds to the identity. 

The aim is to find symmetrized states which diagonal- 
ize the CSCO. It is then easier to combine all the CSCO 
operators into one single equivalent CSCO operator C by 
forming 



C = a H H + a 3 C 3 + a a C a + a s -3(a v i) 



(V.9) 



where the set of factors a must be chosen so that the 
eigenvalues of C are non-degenerate. In the idealized 
case in which the states constituting the solution space 
arc perfect cigenstates of C, non-degenerate a's would 
be sufficient. In practice one should however be careful 
when choosing the factors; guidelines for a good choice 
will be given in Fig. 6. 

Let us now summarize the PTCO procedure. From 
the raw eigenvectors {^^j = l-.-rfs}, a matrix represen- 
tation is constructed for C using Eq. (V.l) and Eq. (V.9). 
The contributing matrix representation for the Hamilto- 
nian in this raw basis is just a diagonal matrix with the 
corresponding raw eigenvalues. The unitary matrix V, 
consisting of eigenvectors of C will enable the diagonal- 
ization towards the diagonal matrix A 



is nothing else than the new corrected energies. One 
notes the presence of the index fi: the PTCO cannot 
restore completely the degeneracies lifted by the original 
numerical implementation. Nevertheless we still expect 
closely packed eigenvalues for each T, I subspace. Second, 
the average values of the C3V class operators are denoted 



(V.13) 



which are closely packed around the ideal value Af 
Eq. (V.7). Again there is a slight fj, dependence. By 
comparison with Eq. (V.7) one obtains the characters \i 
for every state, from which, via the character table, one 
can deduce the irrep T. Third, the corresponding value of 
the canonical subgroup class operator d{a v \) will provide 
the label /1. 

It is trivial to automate also the symmetry recogni- 
tion step. Alternatively it can be performed using the 
projection operators Eq. (V.2), with the advantage of 
giving a number which can be interpreted as a degree 
of symmetry, which is a figure of merit of the achieved 
symmctrization. 



C. A hierarchy of symmetries 



V ] CV = K. (V.10) 

The matrix elements of V provides in one automated 
run the new symmetrized eigenvectors (expressed in the 



With the procedure we have developed, the numerical 
results will first be symmetrized w.r.t. C3V, construct- 
ing the CSCO given explicitly in the previous section. 
The C3 V symmetry is the true symmetry of the structure, 
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only broken by the unsymmctric grid, and numerical er- 
rors. From the physical point of view, there might also be 
higher approximate symmetries like Cg v , D3h and Dgh. 

The C6v symmetry group is constructed from C3V by 
adding mirror operations w.r.t. three intermediate ver- 
tical symmetry planes (<Jdi,i = 1 . . .3), and D3h is con- 
structed from C3V by inclusion of the mirror operation 
of the horizontal symmetry plane (o^). Combining the 
elements of C6 V and leads to the mesoscopic sym- 
metry Dgh. In both Cgv and D3h, there are three double 
group irreps, E1/2, E3/2 and E5/2, all two-dimensional. Bandgap 3 
Their identical subduction tables to C3 V yields the cor- 



respondence E; 



3/2 



E 



3/2 



l E % 



^3/2 1 E b / 2 E1/2 and 
Ei 1 2 — > Ei J 2- The subduction is similar for the Dgh 
group in which there arc 6 similar double group irreps 
of gerade and ungerade type. The irrep matrices D T ex- 
panding the set given in Ref. 1 are given for all groups in 
Fig. 19. 

We may investigate the validity of approximate sym- 
metries by attempting to diagonalizc a new CSCO con- 
structed with respect to the elevated groups. In such a 
step, Cg v , D3h or Dgh class operators are included in the 
set of operators to be diagonalized. This implies choosing 
operators for the Cg = {Cq , Cg~} or the S3 = {S^ , S^} 
classes for Cg v and D3h, respectively 26 . These operators 
enable distinction between Ei/ 2 and £5/2 irreps. For Dgh 
we include both of them simultaneously. 

One thus sees the versatility of the PTCO, which al- 
lows not only to restore the true symmetry, but also to 
investigate a possible hierarchy of symmetries as ever 
cruder, but informative, approximations to the problem 
at hand. 



VI. NUMERICAL RESULTS 

We consider a GaAs QD as described in Sec. II. The 
length of the hexagon edges is 20 nm and the axial length 
of the dot is 5 nm. An Alo.3Gao.7As shell of thickness 10 
nm surrounds the QD in the transverse direction, with 
infinite potential outside. In the axial direction, the dot 
is surrounded by a thick layer of Alo.3Gao.7As, sufficient 
to ensure that all probability distributions goes to zero. 
The parameters utilized in the k-p Hamiltonian are sum- 
marized in Tab. I. 



Numerical implementation, convergence and 
accuracy 



TABLE I. Parameters for k ■ p simulation, at OK. 



Effective el. mass 



Luttinger 



GaAs: m* = 0.067 
AlAs: m* = 0.15 

71 72 73 
GaAs: 6.9800 2.0600 2.9300 
AlAs: 3.7600 0.8200 1.4200 



E GaAs = 1 5lg e y 
E MAs = 3 099 e y 

jpAl x Gai_ :E As 

= (1 - x)Ef aAs + xE, 
c = -0.127 + 1.310a; 



x(l 



Ratio: 60.4:39.6 



34 



Band offset 



E v = 0.396£V 



E c = 0.604S, 



AlxGa±^ x As 

g 

Al x Gai_ x As 



a The Luttinger parameters for the Alo.3Gao.7As is found using a 
linear interpolation of the valence band effective masses and the 
degree of anisotropy, as proposed by Vurgaftman et ai. 33 



which is the energy deviation relative to values obtained 
at highest resolution: 



2m Q m* Vba 



(VI.1) 



Vbar is the potential difference between the dot and the 
first barrier, and L cc i is a characteristic length of the com- 
putational domain, defined as the third root of the com- 
putational volume. For the valence band, m* is replaced 
by a directionally averaged HH mass (Eq. (III. 6)). Using 
A num to characterize the convergence allows to extrapo- 
late by comparing the resolution dependence for the con- 
duction and valence band. 

Figure 4 shows A num for valence band levels 1 and 
15, and for conduction band levels 1 and 17 which con- 
tains similar dominant envelope functions. Only some 
low-resolution points are included for the computation- 
ally demanding valence band, confirming the convergence 
dependence. For a resolution of 91 points, we have 
A num = 1.25 • 10~ 8 for conduction band level 17. This 
corresponds to A# = l.lmeV, and the corresponding 
extrapolated value for the valence band is -0.3 meV. Fig- 
ure 4 also demonstrates that the relative energy differ- 
ences within each band are more accurate than the ab- 
solute values. 



The potential is defined on a square grid using 91 
points in the x, y and z direction, leading to the cross- 
section shown in Fig. 3. On the lateral sides wee see 
the deviation of the numerical implementation with re- 
spect to Cgv symmetry. Any effect from this symmetry 
breaking will later be ideally averaged by the PTCO. 

To assess numerical convergence, we define the dimen- 
sionlcss numerical deviation A num . It is related to 



B. Eigenenergies 

The potential profile of the QD is shown in Fig. 5a). 
In Fig. 5b) we show the calculated lowest order energy 
levels for both electrons and holes. All energies are with 
respect to the top of the valence band. Note that the 
highest valence band level included here deviates from 
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Potential at z=0 



Potential [meV] 



Energy [meV] 




FIG. 3. Cross section of the CB and VB potential in the 
transverse plane. Note the effect of the square grid on the 
lateral edges. 
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FIG. 5. The potential profile along the axis of the nanowire 
(a), and the 10 or 15 lowest order energy levels of the conduc- 
tion and valence band, respectively (b). All levels are twice 
degenerate due to time reversal. In Fig. 5 b), the CB levels 
marked with a dotted line correspond to E states and are 
thus four-fold degenerate. Note that the energy scales are not 
equal. 



ing on whether they are even (gerade) or odd (ungerade) 
with respect to inversion. The gerade single group irreps 
are the ID irreps Ai y9 , A2, g , Bi^ g , -82,9, and the 2D irreps 
Ei t9 and i?2,g, similarly for the ungerade irreps(g — >• u). 
A complete summary of energies and symmetry labels is 
given later, in Tab. VI. 



D. Valence band eigenstates 



FIG. 4. Numerical deviation A num (dimensionless) for con- 
duction band levels number 1 and 17, and for the valence 
band levels 1 and 15. 

the top of the valence band by 39meV. In comparison, 
the split-off band is located a distance A so = 341mcV 
from the valence band top. 



C. Conduction band eigenstates 



In the analysis of the valence band states, we will uti- 
lize extensively the PTCO outlined in Eq. (IV. 5). First 
one needs to symmetrize the eigenstates w.r.t. the true 
C3v symmetry of the original problem, as the small lack 
of symmetry due to the grid may strongly influence the 
choice of eigenstates within degenerate subspaccs. The 
PTCO will also give quantitative information about the 
effect of the symmetrization process. In a second step 
we investigate approximate elevated symmetries with the 
same technique. 



Probability distributions \tp n \ 2 for the 10 lowest levels 
in the conduction band are shown in Fig. 6. The full 
spinorial cigenfunctions ip n can be constructed from the 
scalar envelope functions ip n 2A 

Since the states of the conduction band have the full 
D 6 h symmetry in the present model, we label every en- 
velope function in Fig. 6 by Dg^ single group irreps. The 
Dgh irreps are labeled with subscripts g and u, depend- 



1. Analysis of true Cs v symmetry 

The combined CSCO operator, C, used for postsym- 
metrization according to C3V is given by Eq. (V.9); the a 
parameters were left free in Eq. (V.5). From the numeri- 
cal point of view, however, they must be chosen carefully 
to ensure well separated eigenvalues. To this end let us 
consider the eigenvalues of C. It is easy to show from the 
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FIG. 6. Isosurfaces of probability distributions, |V>cb| 2 , for the conduction band states in order of increasing energy. Each state 
is twice degenerate due to spin. The barrier potential of the surrounding shell is V^hcii- The D6h single group symmetry labels 
are specified for each state. A detailed analysis, summarizing energy and symmetry properties, is given later in Tab. VI 



eigenvalue equation 



that one has 



r 



(VI.2) 
(VI.3) 



i— 3.a.s 



where E 



1 1. . i 



and 



arc defined by Eqs. (V.12) 

and (V.13). The ; eigenvalues can be chosen freely, 

as they are directly related to the at coefficients. Since 

we know the ideal analytical eigenvalues Af of the class 

r 

operators Eq. (V.7), and that the corrected energies E ^ l 
will be close to the original uncorrected energies it be- 
comes possible to choose the a coefficients a priori in 
such a way that the CSCO eigenvalues ; will be well 
separated. To this end we consider the approximate spec- 
trum of Eq. (VI.3): 



E 

i— 3,<r 



Oi\ 4 + a s X s 



(VI.4) 



Note that we singled out the last parameter a s which 
is related to the class operator of the subgroup C s dis- 
tinguishing different partner function indices fi: to avoid 
confusion between Af of C3 V and C s we denoted the latter 
eigenvalue X s ,u,- 



Since one of the eigenvalues may be chosen arbitrarily, 
it is convenient to shift C by ohEq, where Eq is the 
ground state energy, and to make z dimensionless by 



r 

setting a H = 1/ max r , A1 ,;(i? M z 
the set of eigenvalues 



Eq). One then obtains 



r E in.i ~ E o 

c n,i - =r 

max E ,, i — En 



(VI.5) 



such that the range of the energy splitting is normalized 
to unity. Here 



r 

T M = a 



(VI.6) 



i— 3,ct 



The explicit choice of is given in Tab. II, the a's 
can be easily obtained by inverting Eq. (VI.6) 

Several considerations have been taken into account in 
our choice of t£. In particular their relative magnitude 
decides which aspects are given highest priority. E.g. 
choosing 3> 1 would emphasize very large separation 
w.r.t. group symmetry labels and enforce symmetry at 
the price of eventually remixing completely the Hamil- 
tonian eigenstates. On the opposite, tJ«1 would give 
negligible weight to symmetry considerations. We have 
chosen an intermediate set of values, where the relative 
step between neighboring energy levels are a bit smaller 
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-3.5 



-1.5 



0.5 



2.5 



C 3v symmetrization 



TABLE II. Parameters r M chosen to separate between C3 V 
irreps and partner functions 



than the difference in eigenvalue for two functions belong- 
ing to different irreps. This choice was found appropri- 
ate for our structure, where the symmetry of degenerate 
states was slightly broken by a non-symmetric grid. 

Having now determined the set a^, we computed the 
C operator in a raw basis comprising the first 30 valence 
band states (including Kramers degeneracy) ordered by 
energy, and diagonalized it. To visualize the effect of mix- 
ing due to eigenstate symmetrization we plot in Fig. 7 the 
norm of the matrix elements of the unitary transforma- 
tion matrix (Eq. (V.ll)). Figure 8 display the change 
in energy AE after symmetrization, confirming that the 
changes are modest. 
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FIG. 8. Changes in energies due to the symmetrization w.r.t. 
C3v The raw energy spectrum is shown in blue, and the 
difference between the symmetrized and the raw spectrum is 
shown in red. 
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FIG. 7. Norms of the matrix elements of the unitary trans- 
formation symmetrizing the eigenstates w.r.t. C3 V . 

At large, the transformation matrix remixes states only 
within degenerate 2x2 subspaccs corresponding to dis- 
tinct Kramers doublets. There are no intertwined blocks, 
but the procedure would work even in this case, sorting 
correctly the partner functions and symmetrize them. In 
the raw states there were also some mixing between cer- 
tain neighboring energy levels, as evidenced by the se- 
quence of energies shown in Fig. 8. In this finer level one 
sees more subtle mixing extending over different Kramers 
doublets e.g. levels 5 and 6. This is a result of the non- 
symmetric grid, and would be cumbersome to correct us- 
ing a stepwise projection operator procedure rather than 
the automated PTCO. Such mixing increase for higher 
excited states, as expected. These include higher spa- 
tial frequencies, and will accordingly be more strongly 
influenced by limited grid resolution. 

After symmetrization, we also analyze the C3 V irreps 
associated with different spinors. Fig. 9. shows the 



double group weight for the 30 lowest order states, di- 
vided again into 15 Kramers degenerate doublets. The 
eigenstates arc very purely symmetrized w.r.t. C3 V 

(w^ (i^j ~ 1), confirming the symmetry of the QD, and 
the postsymmctrization procedure. The small deviation 
for the highest order states is again due to grid inaccura- 
cies. From Figs. 7 to 9, we see that all the states that were 
most strongly symmetry mixed in the raw eigenspace did 
correspond to Kramers doublets with very small energy 
separation. However, there were also other off-diagonal 
mixing leading to the largest changes in energy. 

We proceed to the analysis w.r.t. UREFs and single 
group irreps, which is shown in Fig. 10. The weights 



(tpi), plotted for each HSFB component, tpi 



i = 1...4 



must be distinguished from the UREF weight w b ' a (ip ) 

defined in Eq. (IV.5). Using Eqs. (IV.5) and (V.5) and 
Eq. (48) of Ref. 1, it is easy to see that they differ only 
by a Clebsch-Gordon coefficient: 



ft/3 



c 



r,r t * ; r 



w 



(VI.7) 



We use w T ^{ipi) to be able to check the content of ev- 
ery HSBF component, in particular when it vanishes or 
when the weight is supposed to be equally distributed 
onto separate components. uF.(ipi) has easy interpreta- 
tion, yielding directly the fractions of one given spinor 
component transforming according to the single group 
irreps. For the C3 V case, the two definitions differ by a 
factor v2 for the weight in the central components of 
the E 1 / 2 spinors, but the definitions are equivalent for 
the remaining UREFs. 

Considering Fig. 10, we first note that the UREFs are 
in accordance with Eq. (IV.3). However, there is a strik- 
ing asymmetry between weights of partner functions 1 
and 2 in the \ l E 3 / 2 ) components of the E 1 / 2 states, most 
notably for levels 2 and 3 and 8 and 9. Due to time 
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reversal symmetry, the two | 1 i?3/2) components and the 
two \ 2 E^/ 2 ) components, should correspond to partner 
functions of E with the same weights. Such asymmetry 
can occur due to mixing of states from nearby energy 
levels if the states bear the same symmetry label. In 
this case, they are mixed in a way that cannot be suffi- 
ciently corrected by the C3 V symmetrization procedure. 
This explanation is supported by the fact that neighbor- 
ing E1/2 levels have opposite weight imbalance in Fig. 10, 
and also by their nearby energies in Fig. 8. We interpret 
such mixing as being a result of the imperfect grid and 
quasi-degeneracies. However, the presence of additional 
quasi-degeneracies should alert us of the possible exis- 
tence of an approximate elevated symmetries, possibly 
able to correctly remix the relevant states. 



C 3v irreps. C 3v symmetrized basis 
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FIG. 9. Double group weights w\ 
the C3v group after C3 V postsymmetrization 
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with r, n being single 



group irreps, for each spinor component in the HSBF basis 
after C3V postsymmetrization. 



2. Analysis of elevated Cf,v and D$h symmetries 



The quasi-degeneracy of the first two E 1 / 2 levels hinted 
at the possible occurrence of approximate elevated sym- 
metries. Such symmetries have previously been found 
both in theoretical 24 and experimental 31 ' 35 work on [111] 
oriented C3V heterostructures. Because the mesoscopic 
structure has D6h symmetry, we shall test the possible 
existence of Cgv, E>3h and Dgh elevated symmetries to 
see if there is a hierarchy and, if so, how it is organized. 

Fig. 11 show the weights of the C3 V symmetrized 
states according to irreps of the Cgv symmetry group. 
Very similar results were obtained for D 3 h- We see no 
change for z E 3 / 2 states which are also well defined us- 
ing these elevated symmetry groups. This is in agree- 
ment with the subduction table for Cg v to C3 V which 
yields E Z j 2 — >Ey 2 ®Ey 2 , the same applies for D 3 h- For 
the remaining irreps, the corresponding subduction yield 
Es/2^Ei/2 and Ei/ 2 ^E 5 / 2 . This is in accordance to the 
results of Figs. 9 and 11. However, we see clearly some ir- 
rep mixing which prevents unique labeling of these states 
in view of the elevated symmetry groups Cg v and D3h. 
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FIG. 11. Double group weights w.r.t. the irreps of the el- 
evated C6v group. The eigenstates have been symmetrized 
w.r.t. the C3V group. 
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-^3/2 -^3/2 E l/2 E l/2 E 5/2 E 5/2 
1 ^2 

-3.5 -1.5 0.65 2.65 0.35 2.35 

TABLE III. Parameters for symmetrization w.r.t. elevated 
symmetry groups C3 V and D 3h . 



For closer investigation we shall perform a new sym- 
metrization, as proposed in Eq. (IV. 5). To this end, one 
extends the set of operators in the CSCO with class op- 
erators for classes of Cg v or D 3 h- This can be easily 
done by adding the operators of Cg or 5 3 , respectively. 
There is no need to add further class operators for canon- 
ical subgroups here. This step allows to distinguish be- 
tween the E1/2 and -E5/2 irreps of the elevated groups. 
Of course one cannot expect the output of these elevated 
symmetrizations to be exact since they are not true sym- 
metries of the QD. Our choice of symmetrization param- 
eters is given in Tab. III. 

The parameters were again chosen to separate clearly 
the different irreps. Note however that the separation be- 
tween for irreps that are indistinguishable in C3V has 
been chosen smaller than the separation between those 
corresponding to different C3V irreps, and also smaller 
than 1, giving the Hamiltonian priority over the approx- 
imate symmetry. 

We refer to the new eigenstates as Cg v - and D 3 h- sym- 
metrized, and denote the unitary transformation matri- 
ces from C3V to C6v by V2, and from Cg v to D 3 h by V3. 
The norms of their matrix elements are shown in Fig. 12. 

The transition from C 3v to C6 V (Fig. 12 a)) leaves 
the l E%/2 states practically unaltered, as expected from 
Fig. 11. These states were already well symmetrized by 
the basic C 3v symmetrization step. As expected we also 
see that quasi-degenerate subspaces that had -E1/2 sym- 
metry in C 3v are now remixed, namely levels 2 and 3, 
4 and 5 and also 8 and 9. Furthermore we see in the 
lower right corner of Fig. 12 a), three non-diagonal blocks 
mixing levels 11, 12 and 15 non-trivially together. The 
asymmetric structure of the intermixing blocks might in- 
dicate that a level above level 15 should also have been 
included in the solution space. Block structures like this 
are useful in revealing nearby states of similar symmetry. 

Interestingly, symmetrization w.r.t. C6 V and D 3 h have 
practically the same effect on the eigenstates. This is 
clearly seen in Fig. 12 b), showing the structure of the 
unitary transformation, V3, between these two eigenstate 
bases. The main effect is merely an irrelevant reordering 
of some neighboring states. For this reason we shall not 
any more comment D 3 h separately. Indeed the similarity 
of Cgv and D 3 h indicate that Dgh symmetry will be the 
most relevant elevated symmetry. 

It has still not been demonstrated that the elevated 
symmetry Cgv is in fact a good approximation to the 
true solution of the original Schrodinger equation. In 
Fig. 14, we show the variation of the energy levels due to 
the imposed symmetry Cgv One sees that the changes 
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FIG. 12. Norms of the matrix elements of V2 and V3, con- 
necting the C6v-symmetrized and D3h-symmetrized states to 
C3V- and C6v- symmetrized bases respectively. 



in energy remain very small (AE < 0.6 meV) compared 
to the confinement energies, confirming the continued va- 
lidity of the Schrodinger equation. Not surprisingly, the 
alterations brought on by the elevated Cgv symmetriza- 
tion does however lead to AE an order of magnitude 
bigger than for the C3 V symmetrization (Fig. 8). Clearly 
this is a manifestation that symmetries that are not true 
symmetries of the physical system are imposed. The Cg v 
symmetrization tend to bring closer energy levels with 
intermixed states. We also see that AE is biggest when 
non-degenerate subspaces are remixed, as expected. 

The weights of the C 6v double group irreps for C 6v - 
symmetrized eigenstates are shown in Fig. 13. The 
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C gv irreps. C gv symmetrized basis 
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FIG. 13. Double group weights w^tp ) of the irreps of the 
C6v group after C6 V symmetrization of the basis. 



gle group decomposition of the components of the l E 3 / 2 
states obviously remains the same as in Fig. 10. It is 
more interesting to consider the E 1 / 2 states, especially 
those corresponding to levels 2 and 3, and also the states 
of levels 8 and 9. For these there is a meaningful differ- 
ence. The weight imbalance present in Fig. 10 has essen- 
tially disappeared in the new eigenstate basis (Fig. 15). 
This is a proof that there was indeed a mixing between 
the quasi-degenerate levels due to grid imperfections that 
could not be retrieved by the C3 V symmetrization, since 
the involved states bared the same labels! In the elevated 
Ce v group, since mixed states bear different irreps they 
are easily disentangled. The restoration of the correct 
balance between the weights of partner UREFs bearing 
C3V single group labels is also a clear signature of this 
fact. 
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FIG. 14. Changes in energies due to the symmetrization w.r.t. 
C6v (very similar results hold for Deh symmetrization). The 
raw energy spectrum is shown in blue, and the difference be- 
tween the symmetrized and the raw spectrum is shown in red. 



identification of the Cg v double group irreps subduces 
correctly towards the previous identification with C3 V 
(Fig. 9). The most interesting aspect of Fig. 11 is cer- 
tainly that the first two levels which were characterized 
by E1/2 C3V irreps (levels 2 and 3) are now distinct Cg v 
irreps, with nearly pure character, similarly for levels 4 
and 5. As the two latter levels were not as much degen- 
erate, their energy change is much larger (Fig. 8). The 
same comment also applies to levels 8 and 9, although 
we see that these levels do not seem to be well-described 
by the elevated symmetry. This is not surprising since 
states of higher energy levels are generally more sensi- 
tive to symmetry breaking contributions due to a micro- 
scopic structure. For the lower levels, the remixed states 
are well defined by £q/2 or E 5 / 2 irrep labels, accordingly 
they approximately obey the symmetry of the elevated 
symmetry groups. 

The symmetries of the spinor components of the C6 V 
symmetrized eigenstates are given in Fig. 15 (using only 
C3V irreps, for easier comparison with Fig. 10). The sin- 
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E 3/2,i E 3/2,i E l/2,i E l/2,i E 5/2,i E 5/2,i 

i = g -3.6 -1.6 0.75 2.75 0.45 2.45 

i = u -3.4 -1.4 0.55 2.55 0.25 2.25 

TABLE IV. Parameters for symmetrization w.r.t. elevated 
symmetry group Dsh- 



3. Dgh as the ultimate elevated symmetry ? 

The results from the symmetrization analysis in view 
of the elevated symmetries D3h and C6v revealed a very 
close relation between these symmetry elevations. This 
is a clear indication that the relevant elevated approxi- 
mate symmetry group may in fact be Dgh which collects 
all the symmetry elements of D3h and Cg v - In fact this 
is not entirely surprising as the heterostructure symme- 
try is Dgh, only broken at the microscopic level by the 
crystal symmetry which in principle does lower it to C3 V . 
This motivates a symmetrization w.r.t. the Dgh sym- 
metry, i.e. with inclusion of class operators for both Cg 
and S3 in the CSCO. The double group irreps of Dgh 
carry the same types of labels as C6 V and D3h , however 
there are twice as many as one has gerade and ungerade 
kinds. The parameters for the Dgh irreps are given 
in Tab. IV, the values for the gerade/ungerade were ob- 
tained by adding/subtracting 0.1 to the corresponding 
values in Tab. III. 

In Fig. 16 we first display the norm of the matrix el- 
ements of the unitary transformation bringing the Cg v - 
symmetrized basis to the Dgh symmetrized basis. Ex- 
cept for levels 8, 9, 11, 12 and 15 there are little changes 
apart from trivial reordering of states. The remixing of 
the highest excited states cannot be analyzed due to the 
stringent limitation in the number of states kept in the 
solution space, but the remixing of levels 8 and 9 can be 
clearly interpreted as a further attempt to give distinct 
irreps E 3 / 2 . u and E 5 / 2 ,u to these levels, that were not 
completely separated in the Cg v symmetrization (-E1/2 
and E1/2 in Fig. 11). We clearly see that the depar- 
ture from Dgh symmetry, due to the underlying symme- 
try of the Zinc Blende lattice, starts to be important after 
level 10, in particular for levels 11,12 and 15 which have 
strongly mixed character with respect to Dgh symmetry. 

The energy changes AE due to the Dgh symmetriza- 
tion are given in Tab. V. From this we see again that 
neighboring levels bearing E 1 ^ 2 labels in view of C3V 
are renormalized. The energy differences are very sim- 
ilar to the corresponding values for Cg v (Fig. 14). In 
Fig. 17a), we show the double group weights of the Dgh- 
symmetrized eigenstates. The purities of the states w.r.t. 
the double group irreps displayed in Fig. 17a) are also 
very similar to the corresponding Cg v analysis (Fig. 11). 
The Dgh symmetry is thus fulfilled to approximately the 
same degree as the intermediate symmetries Cg v and D3h, 
as is also confirmed by the deviation in energy due to the 
Dgh symmetrization (Tab. V), which stays on the same 
order as for Cg v . 
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FIG. 16. Unitary transformation between bases symmetrized 
w.r.t. C6v and symmetries. 



The single group analysis of each spinorial component 
in the HSBF basis in Fig. 17 b) leads to results in fair 
agreement with the UREFs for Dgh which are listed in 
Eq. (B.6). These results fully confirm the relevance of the 
approximate symmetry elevation to Dgh for the valence 
band in our structure, and clearly shows when it applies 
and when it has limitations. In the following section a 
more in depth analysis of the UREF weights of Fig. 17 is 
carried out. 



E. Characterization of DPGPS, UREF's and 
envelope function symmetry of valence band 
eigenstates 

Up to now we have concentrated on the symmetry 
properties of eigenstates as a whole. It is however possi- 
ble to get important information and insight within each 
of them by studying the nature of the UREFs, identi- 
fying dominant HSBF (DPGPS) and dominating UREF. 
Such an approach will not only provide enhanced physical 
interpretation and intuition, but can also reveal to be ex- 
tremely helpful information when building more complex 
objects, e.g. multiexcitons in a configuration interaction 
approach and in the strong confinement limit (c.f. 31,35). 

The analysis will be carried out using essentially ele- 
vated Dgh symmetry which provides the finest classifica- 
tion, establishing the DPGPS in Fig. 17 and determining 
the dominant UREF's in Fig. 17. The analysis is finally 
completed in Fig. 17 by an analysis in terms of more clas- 
sical azimuthal and radial quantum numbers, which pro- 
vides an intuitive understanding for the sequence of irreps 
and outlines a link to symmetry adapted functions 26 . 

The main results of this section are summarized in 
Tab. V. Together with Fig. 17, Tab. V provides a deep 
understanding of all valence band states which will al- 
low to decipher the remarkable properties of the optical 
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FIG. 17. Double group weights and single group weights for 
each HSBF component w.r.t. D6h group, using the D6h sym- 
metrized basis. 



transitions. 



1. Heavy and Light hole mixing, HSBFs and DP GPS 

The most common way to analyze valence band states 
in heterostructures is in terms of valence band mixing be- 
tween HH and LH states. This is an import from quan- 
tum well physics where, at zone-center, there is no mix- 
ing. In strongly oblate (c.f. Fig. 1) disk-like structures 
this is still a good starting point due to weak mixing for 
the ground states, allowing to deduce that the ground 
state should be HH-like. Moreover one can predict fur- 
ther the order of excited HH- and LH-like states coarsely 
using a scalar anisotropic effective mass model based on 
Eqs. III. 6 which neglects valence band mixing and as- 



sumes an infinite cylindrical potential well with the same 
height and cross sectional area as the nanowire QD. One 
obtains a separation of about 45 meV between the ground 
HH and LH states. Accordingly, we do expect the 15 low- 
est energy levels (Fig. 5) to be dominated by HH states, 
in agreement with what can be seen in Figs. 10 and 17, 
where it is easy to see that the HH-weight always re- 
mains above 0.5 using Eqs. (IH-7) and (III. 8). We also 
see clearly in Figs. 10 and 17 that the degree of band- 
mixing increases with increasing excitation levels, as ex- 
pected. It is very low for the ground state, 99% HH, and 
starts to be significant above level 6. 

We now proceed to identify the DPGPS, i.e. the dom- 
inating HSBF of a mixed state. Being HH-like, the 
DPGPS of all the first 15 Kramers conjugate pairs is 

#|/ 2 y ,i = 1,2 in C3 V , and the corresponding DPGPS 

in D 6 h are the partners of (#3/2.9) (c.f. Eq. (A. 6)). 

The concepts of DPGPS and dominating UREF are 
important and will enable understanding of the quasi- 
degeneracies in the next section. 



2. Analysis of UREFs 

Let us first identify for every level (Kramers' pair) the 
dominant UREFs, using Fig. 17. We find the following 
dominances: level 1, 6 and 14 are Ai, g , levels 2, 3, 8, 9 
and 13 are E\ tU , levels 4, 5, 11, 12 and 15 are #2,9, level 
7 is i?i. u and level 10 is #2.u- This is summarized in 
column 6 of Tab. V, column 4 (for C3 V ) can be obtained 
by subduction. Note that we have already regrouped 
some levels together as pairs in this enumeration, the 
reason will become clear during the forthcoming analysis. 

The concept of product states between the dominant 
UREF and DPGPS allows to explain most of the level 
sequence in Fig. 17. Indeed, the double group irrep of 
each level (or the irrep pair in case of a level pair) can be 
faithfully generated by making the following products of 
irreps: for the isolated levels 1, 6 and 14 A\ t9 x #3/2,9 = 
#3/2,9, for the level pairs 2, 3 and 8, 9 #i )1( x #3/2,9 — 
E5/2.U + #i/2,mj f° r level pairs 4, 5 and 12, 15 #2,9 x 



E- 



3/2,9 



E, 



5/2,9" 



#1/2,9, for level 7 Bi^x £: 3/2 



for level 10 B 2 u x E- 



3/2,g 



E- 



3/2,9 



#3/2,«, for level 11 E 2 ,, 



3/2, tl, 



#3/2, g = #1/2,9 (+#5/2,9 but the latter is missing), and 
finally for level 13 Ei tU xE 3 / 2g = #5/2, u (+#1/2, u but the 
latter is missing). The missing levels related to level 11 
and to level 13 likely lie above the 15th energy level. The 
association of level 12 (instead of 11) with level 15 will be 
explained in the next section when looking at azimuthal 
and radial quantum numbers. It is still an open question 
whether the departure from Dgh symmetry of levels 11, 
12 and 15 may be further minimized by including the 
missing partner of 11 in the symmetrization procedure 
(c.f. Fig. 16). 

The clustering in energy of the level pairs 2, 3 and 4, 
5 and 8, 9, and 12, 15 (see Fig. 14 and Tab. V) can also 
be very clearly explained by dominant product states. 
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It should be kept in mind that the quasi-degeneracies 
of 4,5 and 12,15 with levels 6 and 13,14 respectively are 
accidental since these have different symmetries. Whilst 
the energy of the lower pairs are closely packed together 
due to negligible valence band mixing, higher pairs like 
8, 9 and 12, 15 become more and more significantly split 
apart. The splitting is however already visible for the 
level pair 4, 5, manifested in Fig. 17 as non-vanishing 
weight in the central E x i 2 components in which there are 
distinct UREFs for the respective levels. This is a clear 
demonstration that band-mixing is also a very relevant 
concept in the HSBF picture. 

The isolated levels 1, 6 and 14 are all linked with the 
dominating A\ t9 UREF. Although it is customary that 
the fundamental level is strongly dominated by a fully 
invariant envelope function concentrated in the DPGPS 
component, the structure of levels 6 and 14 can again 
only be explained in the next section when looking at az- 
imuthal and radial quantum numbers. Note that level 14 
also display contributions from other HSBF components. 

Levels 7 and 10 are also interesting. They are very 
similar and in a way form a complementary pair because 
the relative strength of the weights of B\ u and B 2u is 
opposite between the two levels. In fact within one level 
none of the two is strongly dominant, but this is allowed 
by symmetry since the two products B\ M x E 3 / 2g and 
B2, u x E 3 / 2 .g gi ye both of them global E 3 / 2 ,u symme- 
try. In contrast with the partners of an E lu / g or E 2u / g 
irrep (see Fig. 17) which must be always balanced as 
predicted by the Wigner-Eckart theorem, the weights of 
the Bi u and B 2u UREFs must not necessarily be bal- 
anced in these levels. Levels 7 and 10 therefore illustrate 
the general fact that even if a state has a clear DPGPS 
(clearly dominant HSBF) it may not necessarily have a 
single very dominant UREF, as allowed by the analytical 
UREF decompositions given in Eq. (B.6). 

Thus far we have only considered the DPGPS compo- 
nents of the spinors, however the higher excited levels 
also have a significant contribution from the E 1 i 2 HSBFs 
(LH) . For simplicity we shall discuss this only in C3 V sym- 
metry. First, note that all i?-partner functions within LH 
components have balanced (equal) weight despite band 
mixing: this is imposed by symmetry (c.f. Eq. (IV. 3)). 
By contrast, and in analogy to the discussion of the previ- 
ous paragraph, level 13 has A i} i — 1,2 UREF pairs in the 
LH components with an allowed weight imbalance. Sec- 
ond, we have seen analytically that time reversal sym- 
metry imposes restrictions (c.f. Eq. (IV. 4)) on the LH 
components of the £1/2 states. As a clear example let 
us consider the first pair of E 1 / 2 states with significant 
LH weights, levels 8 and 9, which in addition contain 
all UREFs predicted by Eq. (IV. 3), as can be seen in 
Fig. 11. We have been able to check numerically that in- 
deed the relative phase of the A\ and A 2 UREFs are out 
of phase with the E UREFs by a factor i in agreement 
with Eq. (IV. 4). Note also that level 12 has a rather large 
Ai LH component, we do not think however that this 
state should be identified with a ground LH-like state. 



Finally, it is interesting to consider the parity of the 
dominant UREFs found with respect to the operation. 
From the T)eh character table it is easy to see that all 
the single group irreps of the dominant UREFs appear- 
ing so far, namely A lt9 , A 2 , g , B lyU , B 2 , u , E 2 , g and E ltU are 
even under ah (nevertheless it should be noted that non- 
dominant UREFs may acquire some level of z excitation 
due to bandmixing (c.f. Fig. 17)). The reason for this 
behavior of the dominant UREFs is rooted in the high 
lateral to axial aspect ratio of the QD under consider- 
ation which is 40:5 (c.f. Fig. 1). We therefore expect 
the envelope functions of the 15 lower levels considered 
here to have excitations mainly in the transverse plane 
(azimuthally and radially). This is investigated in the 
following subsection (Fig. 17). 

To conclude we would like to stress that the dominant 
UREFs that we have identified in this subsection will 
play a key role not only in enabling the use of azimuthal 
and radial quantum numbers, but also in analyzing the 
strength of optical transitions. 



3. Approximate azimuthal and radial quantum numbers: 
the role of symmetry adapted functions 

The goal of this section is to investigate the correspon- 
dence between a classical effective mass analysis for a 
disk-shaped quantum dot in terms of standard azimuthal 
and radial quantum numbers (as we have just seen there 
is no need for a vertical n z quantum number), and the 
sequence of irreps for the dominating UREFs identified 
in the previous section. 

The Schrodinger equation for a circular disk-shaped 
QD with infinite barriers is separable in azimuthal and 
radial coordinates, and leads to eigcnfunctions which 
are proportional to products of an azimuthal expo- 
nential exp(im a (f)) with an Z a -th order Besscl function 
Ji a (ki a ,n r r ) f° r the radial coordinate, where it is un- 
derstood that ki a . nr R is the n r -th zero of Ji a with R 
being the dot radius. They are labeled by the az- 
imuthal and radial quantum numbers, l a = 0, 1, . . . , 00 
and n r — 1, . . . , +00. In cylindrical coordinates, m a is 
restricted to m a = ±l a (and axial excitations are labeled 
by n z ). The number of azimuthal nodes is thus l a , and 
there is a twofold degeneracy for states with l a ^ 0. 

Similar approximate quantum numbers can be ex- 
pected in our nearly cylindrical 3D solutions, see Fig. 6 
for the conduction band, and Fig. 18 for the valence band 
spinor component cross sections. Therefore, in the fol- 
lowing we shall specify the character of the computed 
dominating UREFs not only by their irreps, but also by 
additional subscripts ipn a n r where n a and (n r — 1) spec- 
ify the number of nodes in azimuthal and radial direction 
respectively (when they can be determined). The sub- 
scripts were identified by visual inspection for the first 
15 levels, and are listed in the last column of Tab. V. 

We now consider the energy level sequences. First the 
sequence of levels 1, 6 and 14 with dominating UREF of 
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FIG. 18. Cross sections of a selection of the contributions 
to the valence band spinor. The envelope functions shown 

A, A, 

in the top row are characterized as ip 01 9 and ip 02 9 , in the 
order from left to right. The two functions in the lower row 

,Si,„,l , ,-Bi. u ,l 

are tp 11 and ip 12 



A\ t g symmetry: it is easy to see that they nicely form the 
sequence of m a = levels with n r = 1,2,3, i.e. purely 
radial excitations. 

Another less trivial, but fundamental, sequence is the 
n r = 1, l a = 1 . . .5 sequence. We recognize the follow- 
ing five level pairs: 2 and 3, 4 and 5, 7 and 10, 11 and 
missing partner, 13 and missing partner. The question is 
now how to obtain further confirmation of this sequence. 
If one looks at the subduced representations from 0(3), 
the so-called Symmetry Adapted Functions (SAFs), listed 
in Table T35.5 and T35.6 of Ref. 26, we can confirm 
that the link is correct, as the SAF quantum numbers 
I, to, here corresponding to Z Q , m a are in agreement with 
the appropriate single group labels. One finds that the 
I = l,m = ±1 basis subduces to E% jU in Dgh symmetry, 
I = 2, to = ±2 subduces to E 2 „, I — 3, to = ±3 subduces 



and 



to i?i jtI + i?2,u, I = 4, to = ±4 subduces to E2, g , 
finally , I = 5, m = ±5 subduces to E\^ u . Hence we have 
confirmed that this sequence of azimuthal excitations is 
in agreement with the sequence of irreps for the dominat- 
ing UREFs, and additionally that levels 11 and 13 have 
missing partners which must be at higher energies. For 
level 11, this was also indicated by the structure of the 
intermixing blocks in the lower right corner of Fig. 17 as 
previously mentioned. 

Similarly we can identify the last n r = 2, l a = 1,2 
sequence corresponding to level pairs 8, 9 and 12, 15. 
Again, the corresponding to quantum numbers subduce 
to the correct irreps E\ tU and E 2 , g respectively. 

Even further interpretation is sometimes permitted. 
Let us consider again the pair of levels 7 and 10. Three 



nodal planes here correspond to either B\ u or B 2u sym- 
metry. These are distinguished by a rotation by 7r/6 of 
the nodal planes. In our choice of coordinate system 
(Fig. 2), the ipsx' u functions have nodal planes perpen- 

lB 2 ,u 



,Bi,« 
1 

dicular to the facets of the hexagon, whereas the ip^i' u , 



have nodal planes intersecting the corners of the hexagon. 
These two orientations give rise to different filling factors. 
The i/j^i'™ functions fill the area of the hexagon more ef- 
fectively than the corresponding functions with B2.U la- 
bels. The curvature, and thus the energy is therefore 
higher for the ipa''" ■ This is also seen in the conduction 
band, as there is a splitting between the two correspond- 
ing levels (CB levels 7 and 8). For the valence band, the 

spinors do in fact contain mixtures of ip3i ' u arL d V , 3i 2 '"- 
Energy level 7 has B ljU as the largest contribution and 
a smaller contribution of B 2 ^ u . The situation in energy 
level 10 is mirrored. The energy splitting between the lev- 
els is mostly due to the different curvature of ip^x' u and 
il>3l' u ■ Note that the different curvature of states rotated 
by 7r/4 also apply for the E1/E2 states with azimuthal 
excitations. However for these states unequal weight in 
the UREF partners would violate the restrictions from 
Wigner-Eckart theorem. 

As valence band mixing is weak for the lowest energy 
levels, the corresponding dominant HH envelopes will not 
hybridize with LH components. The lowest order UREFs 
observed in the LH components are ipfi' u ,i — 1..2 enve- 
lope functions. This also further supports our previous 
assumption regarding the DPGPS as there is no sign of 
an LH-like sequence. 

Note that the considerations of this section can also 
be carried out for the simpler conduction band, we leave 
this as an exercise for the reader. The situation in the 
valence band is much more complicated due to direction- 
dependent effective masses and valence band mixing. We 
expect that the approach also works when sequences of 
vertically excited states occur, and with E 1 / 2l9 DPGPS 
states (LH-like). 

At this stage Table V entails the full information for a 
complete classification of all states in compact form, both 
in true C3V symmetry and in elevated Dgh symmetry. 
Note that the energy deviation due to symmetrization is 
the same in D@h as in C3V for the states labeled E 3 / 2 ,g 
and E 3 / 2 ,u (Dgh)- Also note that in Dgh, the energy 
deviation is almost balanced within one pair. The last 
column of the table summarizes also the azimuthal and 
radial properties of the dominating UREF. 

In Tab. VI we present a similar summary for the con- 
duction band states. 



VII. OPTICAL TRANSITIONS 

Sharp optical transitions in quantum dots are among 
their major attractive features for device applications. It 
is well-known that the oscillator strength corresponding 
to each interband d-polarized transition (d = x, y, z) at 
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TABLE V. Final classification of valence band states. 



Level 


^raw 

[meV] a 


A£ C3 - 


A£ D6h 
[peV] 


«& 1]b 




UREF° 3v 
dominant 


0» D6h 

dominant 


UREF Dfih 
dominant 


' (n a ,n r ) 


1 


-14.6 


1.3 


1.4 


0.01 




Ai 


E 3/2,g 


Al,g 


A i 

V> i 1,3 


2,3 


-20.4, 
-20.9 


0.4, 
0.5 


118.2, 
-111.5 


0.04, 
0.03 


E 1/2 




E$/2,u: 
E\/2,u 


E\ tU 




4, 5 


9fi A 

28.0 


A 8 

14.9 


-487.5 


n 1 9 
0.08 


El/2 




Eb/2, 9 , 
E l/2,g 


E 2,g 




D 


98 K 


-10. D 


-io. y 


U. 10 


E 3/2 




&3/2,g 


A i,g 


^02 


7, 10 


-31.4, 
-35.0 


6.3, 
-9.5 


8.9, 
-12.1 


0.25, 
0.21 


E j 

^3/2 


Ai, 
A 2 


Ez/2,u 


Bi, u , 

B2,u 


V>31 , 

V'si 


8, 9 


-33.9, 
-34.6 


-2.8, 
8.5 


319.8, 
-296.3 


0.32 


-&L/2 




El/2,u: 
Es/2,u 


Ei :U 




11 


-35.4 


1.2 


285.0 


0.37 


E 1/2 


E 


C El/2,g 


E2,g 




12, 15 


-37.8, 
-39.4 


7.5, 
-29.5 


-57.7, 
-281.8 


0.42, 
0.43 


Ei/2 


E 


c Ei/2,g: 
E 5/2,g 


E 2,g 




13 


-38.6 


-3.5 


-27.1 


0.46 


E 1/2 


E 


C E 5 /2,u 


Ei iU 


</> 5 i ' 


14 


-39.0 


16.3 


15.7 


0.44 


E 3/2 


A 1 


E 3/2,g 




V^03 



a Averaged over Kramers doublets 
b For the C3 V symmetrized basis 
c Intermixed symmetries 



TABLE VI. Final classification of the conduction band states. The (r, /x) D(ih labels are obtained taking the direct product of 
the UREFs and the DPGPS representations (recall that spin-splitting is ignored). The DPGPS for all conduction band states 
is S 1/2iU (c.f.Eq. (A.6)) 



Level 


^raw 

[meV] 


(r, M )°3v 


UREF 

C3v 


(r,M) Deh 


UREF 
D 6h 


w, (r)D6h 

^(n a ,n r ) 


1 


1590 


-&L/2 


Ai 


El/2,u 


Al,g 


<i' s 


2,3 b 


1599, 
1599 


^3/2' 
El/2 


E 


E 3/2,g, 
El/2,g 


Ei zU 




4,5 b 


1612, 
1612 


El/2, 

E J 

^3/2 


E 


E$/2,u, 
E3/2,u 


E 2,g 




6 


1617 


El/2 


Ai 


El/2,u 


Al,g 


Wd2 


7, 8 


1626, 
1630 


El/2 


Ai, 
A 2 


E 5/2,g 


Bi,u, 

B2,u 


^31 , 

,B 2 , U 


9, 10 b 


1636, 
1637 


El/2, 

E j 

^3/2 


E 


El/2,g, 
E 3/2,g 


Ei zU 


V>12 



Degenerate by symmetry. To enable better comparison with the valence band, we do however preserve individual level numbering for 
all 10 conduction band levels despite the degeneracy. 
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frequency hu> = E^ c — E^" is proportional to the summed 
squared c — v interband matrix element Afd i „ cj „ v (r c ,r v ), 
where the sum is over the degenerate contributing states, 
and can be written 

M d ,„ c> „ v (r c ,r v )= (vii.i) 
EE E : n ^K^c\ Pd IVv : n v ,r;,^ v )| 2 . 

The primed sum means a sum over eventually conjugated 
contributions, i.e. T' summed over {T, V*} only if V ^ 
r*. In Eq. (VII.l) ip c and ip v denote full conduction and 
valence band kets corresponding to energy levels E^ c and 
respectively. 

For the symmetry groups considered here, the set of 
momentum operators P4, d = x,y, z, can be divided into 
two distinct irreducible tensor operator (ITO) sets. In 
C3V one has 26 : 



{Px,P„}<~ 

and in Dght 

\Pxi Py} ' 



{^f,/i=l,2} , P z ^rP A - (VII.2) 



/i 



1,2} 



P z 



The correspondences (VII.2) 
the use of the generalized 
orem for evaluating the inner 
(■0c : n c ,T c ,/j, c \ P d |0 V : n v ,T v ,/j, v ). 



, pA 2 , u 
(VII.3) 

and (VII.3) allow 
Wigner-Eckart the- 
matrix element 
For our symme- 



the prediction of optical 
in the x — y plane 24 , and 
selection rules for optical 



try groups this leads to 
oscillator strength isotropy 
to use group theoretical 
transitions. 

It should be stressed that the decomposition into 
UREFs linked with the HSBF basis allows to go beyond 
the simple use of the Wigner-Eckart theorem. Indeed, 
the UREFs lead straightforwardly to the prediction of 
"magic ratios" in polarization anisotropy , in the frame 
of our Luttinger model, and a conduction band with no 
spin-splitting. To this end, full conduction band spinors 
must be reconstructed by adding spin as in Ref. 24. In 
the case of C3V symmetry, four types of Kramers degener- 
ate pairs then appear : E 1 / 2 (Ai) and E 1 / 2 (A 2 ) for con- 
duction band states with A\ and A 2 envelopes respec- 
tively, and Ei / 2 (E) and E 3 / 2 (E) f° r conduction band 
states with E envelopes. The latter two are always de- 
generate L Then, thepolarization anisotropy, defined as 
Aij = (Mi - Mj)/(Mi + Mj),i,j e {x,y, z} can be cal- 
culated using UREFs, and reveals the magic anisotropy 
ratios A zx given in Tab. VII. The ratios ±1 in Tab. VII 
stem from standard C3V selection rules. The Wigner- 
Eckhart theorem also predicts A xy = 1 and A zy = A zx 
for all irreps. 

In the present work we shall now show that when the 
identification of the dominant UREF can be performed, 
one obtains further clues to the oscillator strength spec- 
trum. Moreover, when azimuthal and radial quantum 
numbers can also be identified, further constraints occur. 
In both cases the presence of main lines can be explained, 
and fine-structure can be interpreted. 



TABLE VII. Optical transition anisotropy ratios, A z 
QD's with C3V symmetry 



for 



"\ r v 
r c ^\ 


l E 3 /2 Ei 12 


l E 3 / 2 (E) 
Ei /2(E) 
E 1/2 (Ai) 


1 -1 

-1 f{^^,^T 
-1 ! 



Numerical calculation of C3V optical oscillator 
strengths 



The calculated optical oscillator strengths of the 
nanowirc QD arc displayed as a function of the transi- 
tion frequency in Fig. 19. We used D@h symmetrized 
valence band states since it allows to interpret the opti- 
cal spectra in most details. In the side figure the details 
of two transitions calculated in the C3 V -symmetrized ba- 
sis are also displayed, other differences were minor. By 
contrast, for the conduction band we used the raw cal- 
culated states since these states were already sufficiently 
symmetrized, rendering unnecessary a further D 6 h sym- 
metrization. Nevertheless we restored in Fig. 19 an exact 
degeneracy for E conduction band states, for the purpose 
of clarity of the fine structure in the figure. The numeri- 
cal splitting of the computed energies was anyway below 
the estimated convergence (cf. Fig. 4). 

The upper subplot in Fig. 19 shows the total oscilla- 
tor strengths for the main transitions (summed over in 
all directions), whilst the lower subplot shows the corre- 
sponding A zx anisotropy ratios. This data suffices for our 
purpose, but the reader may recover the separate results 
in each direction using the analytical results A xy = 1 and 



A zx , whose validity was also confirmed numeri- 
cally. It must be understood that other effects also oc- 
cur in a real experiment, besides the neglected Coulomb 
contributions. They may change to some extent the pre- 
dicted spectrum, in particular in nanowires there are im- 
portant effects due to the high index contrast between 
the nanowire and the air 36 (dielectric mismatch effect). 
Nevertheless, the oscillator strength spectra given by sin- 
gle particle calculations as in Fig. 19 are often the main 
characteristic of the intrinsic optical response of QDs. 

Let us now study more in details the oscillator 
strengths of each optical transition, which is color coded 
in Fig. 19. They are numbered so that the properties 
of the conduction band level i and valence band level j 
corresponding to transition CBi — VBj can be directly 
read off using the previously obtained Tabs. V and VI, 
and Eqs. (VII.2) or (VII.3).' 
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FIG. 19. The top subplot show the combined optical oscillator strength, M x + M y + M z , as a function of energy. The colors label 
the contributions due to the individual conduction band levels, and the valence band state number is indicated by numbers. The 
lower subplot show the anisotropy ratios, A zx , for the significant transitions (limit set to M x +M y +M Z =0.05). The theoretical 
anisotropy ratio predictions are also included. Some of the anisotropy ratios cannot be predetermined, for these transition the 
theoretical value has been set to the dummy value to mark them. The inset at the right show the CBi — VBj, i,j = 2, 3 in 
the C3v basis. 



B. Dominant and missing optical transitions 

A first striking feature of the upper part of Fig. 19 
is the dominant diagonal character of the optical transi- 
tions, i.e. CBi — VBj is most intense for i = j. Note that 
within degenerate CB levels the numbering has been ex- 
plicitly chosen to respect CBi — VBj diagonality to ease 
the analysis. This diagonal character is nearly perfect 
for the ground state transition, as well as for the first 
pair i = 2, 3, but progressively weakens as one climbs the 
excitation ladder. 

We also clearly see in Fig. 19 that the highest pairs 
give rise to a richer structure with side peaks due to the 
valence-band complexity, but still the (CB7, CB8) pair 
is dominantly diagonally coupled to the (VB7, VB10) 
pair, and (CB9, CB10) to (VB8, VB9) pair (see Tab. V). 
Clearly, cross-coupling like CB7-VB10, e.t.c. also appear 
as sub-structures within these dominant pairs. 



Looking back to the character of the corresponding 
states in Tabs. V and VI, we find that symmetry eleva- 
tion to Dgh is an important ingredient to explain missing 
non-diagonal transitions in Fig. 19. However, not all of 
them are explained. Take for example the set of transi- 
tions from the ground conduction band level CBI — VBj. 
Only CBI — VB1 is observed. In principle, using double 
group selection rules in C3 V , all values of j are permit- 
ted from CBI in x — y-polarization, but in Dgh only 
j = 1,5,6,11,12,14 are permitted. In z-polarization 
j = 2 — 5,8,9,11 — 13,15 are permitted in C3 V whilst 
only j = 5, 11, 12 would be in Dgh- Thus symmetry el- 
evation accounts for part of the missing transitions, but 
obviously it is not enough to explain all of them. 

Actually one should consider in addition both the sym- 
metry of the dominant UREF and the presence of the ap- 
proximate azimuthal and radial quantum numbers to ex- 
plain all missing transitions for the lowest energy states. 
This is particularly clear with the set CBI — VBj. First, 
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on the basis of the symmetry of the dominant UREF 
A 1>g of CB1 and the contributing UREFs of VBj, all 
z-polarized lines are forbidden, and only j = 1,6,14 
are allowed in x — y polarization. Second, to explain 
the remaining missing lines from CB1, corresponding to 
j = 6, 14, it suffices to invoke the approximate selection 
rule based on the radial quantum number of the domi- 
nant UREF, which is n r = 1 for CB1, and n r = 2, 3 for 
VB6 and VBU respectively (see Tabs. V and VI). 

Let us examine more closely the optical coupling be- 
tween diagonally coupled pairs, considering the UREFs. 
In (CB2,CB3) - (VB2,VB3), symmetry elevation to 
D 6h is not sufficient to predict only CB2 — VB2 and 
CB3 — VB3 transitions as observed in Fig. 19. In this 
respect, it is interesting to see that if one assumes C3 V 
symmetry, both the crossed transitions between the pairs 
(CB2-VB3 and CB3 — VB2) should appear, and indeed 
this is what we observe in the calculated spectrum in the 
C3 V -symmetrized basis (see side of Fig. 19). The obser- 
vation of CB3 — VB2 would witness the Dgh-symmetry 
breaking present in the [111] Luttinger Hamiltonian, and 
might well be resolvable experimentally, whilst CB2-VB3 
would be compatible with both symmetries. 

In the higher diagonally coupled pairs (CBA, CB5) — 
(VB4, VB5) and (CB9, CfllO) - (VB8, VB9) the same 
effects are observed and the same considerations ap- 
ply. For these two sets of transitions we also see in 
Fig. 19 the appearance of additional non-diagonal transi- 
tions with other levels, a trend which naturally increases 
with transition energy. The last diagonally coupled pairs 
(CB7, CB8) - (VB7, VB10) are particularly interesting 
because the two electronic levels CB7 and CB8 are non- 
degenerate, hence in Fig. 19 one observes four dominant 
peaks for these transitions. 



C. Optical anisotropy of dominant transitions 

The optical anisotropies are also of fundamental inter- 
est. By contrast to the selection rules discussed in the 
previous section, nearly all optical anisotropies can be 
explained in C3 V symmetry. 

Table VII, with the help of the classification of Tabs. V 
and VI, predicts the anisotropy ratio to be —1 for tran- 
sitions between isolated levels CBi — VBi, i = 1 or 6, in 
agreement with the lower part of Fig. 19. 

The transitions between the coupled pairs (CB7, CB8) 
and (VB7, VB10) are also predicted to have anisotropy 
ratio —1, as is well verified in Fig. 19. However the situa- 
tion for the other diagonally coupled pairs is more subtle. 

In (CB2,CB3) - (VB2,VB3), the optical anisotropy 
is predicted to be undetermined for the two transitions 
stemming from CB3, and —1 for the two stemming from 
CB2. This is satisfied in the C3 V -symmetrized basis (see 
side of Fig. 19), but the anisotropy of both, particularly 
CB3 — VB3, is still close to —1. The explanation is 
very simple if one recalls that the LH weight in VB2 and 
VB3 is less than 4% (from Tab. VII, the undetermined 



anisotropy ratios will be close to -1 if the overlap with 
$^ vanishes). The departure from —1 is much bigger for 
CB3 — VB2, which can be understood by recalling that 
this transition should be forbidden in the x — y polar- 
ization in Dgh , hence it is more sensitive to LH admix- 
ture. We should point out that the disappearance of the 
transition CB2 — VB3 in the D 6 h-symmetrized basis, ev- 
idenced in Fig. 19, cannot be understood on the basis of 
Dgh symmetry alone. We suspect that a study general- 
izing the magic ratios to D 6 h symmetry, linked with the 
neglect of conduction band spin-splitting, might explain 
this observation, but the analytical verification would be 
quite overwhelming. 

The optical anisotropy of (CB4, CB5) - (VB4, VB5) 
and (CB9, CB10) - (VB8, VB9) is entirely similar to 
(CB2, CB3) - (VB2, VB3) but with the role of the two 
conduction and valence band levels reversed. We also 
see that the undetermined anisotropy ratio in transitions 
from CB9 to (VB8, VB9) depart more significantly from 
any specific values, as is expected due to increased band 
mixing. 



D. Fine structure due to valence-band mixing 

Besides the series of dominant peaks there is a fine 
structure produced by valence band mixing which is par- 
ticularly interesting to investigate. Fig. 19 shows that 
such band mixing really starts to be significant from the 
set of optical transitions CBi — VBj with i, j > 4. 

From the pair (CBA, CB5), weak additional transi- 
tions to VB6 and VB11 are visible in Fig. 19. From 
the Wigner-Eckart theorem, the transitions to VB6 are 
both allowed in elevated symmetry D 6 h whilst towards 
VB11 only CBh - VBW is allowed. We conclude that 
CBA — VB11 is a manifestation of the true C3 V dot sym- 
metry. These weak transitions are slightly more intense 
in ^-direction. The anisotropy ratios for the transitions 
CBA—VB6 and CB5-VB11 arc -1, and for CB5-VB6 
it is +1, in agreement with C3 V symmetry (Tab. VII). 
For the remaining CBA — VBW transition the undeter- 
mined anisotropy ratio proves to be numerically close to 
+1. Clearly, looking at the corresponding i?2. 3 UREFs, 
in the spinors Eqs. (C.lc) and (C.ld) we may confirm 
that both +1 transitions are purely due to valence band 
mixing since only the ^1/2. 3 HSBF (LH-components) can 
couple to CBb. 

Besides the next dominant transition CB6 — VB6, we 
can see two weak transitions CB6 — VBA and CBQ — VBh 
in Fig. 19. Again, both transitions are due to valence 
band mixing, and are not predicted in elevated symme- 
try D 6 h; they manifest the true C3V symmetry of the QD. 
The novelty here is that both transitions should have the 
"magic" anisotropy ratio 3/5 predicted by Tab. VII. This 
is very well satisfied by CB6 — VBA, but CB6-VB5 dis- 
plays a numerical anisotropy ratio close to 0.3 which will 
represent the largest deviation from the predicted value 
in our numerical data. Since its oscillator strength is al- 
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ready rather weak we can safely attribute this deviation 
to imperfect C3 V symmetrization. 

Let us now investigate the weak transitions from the 
(CB7,CB8) pair to levels VBj, j = 8,9,13. Again 
these six transitions due to valence band mixing man- 
ifest the true C 3v symmetry of the QD, and display 
quite accurately the "magic" anisotropy ratio 3/5 for 
£-1/2(^1) — -E-i/2 transitions. 

The last weak transitions are related to the highest 
excited pair (CB9, CB10) towards VB7 and display po- 
larization anisotropy ratios of (+1, —1) respectively. The 
symmetry considerations are the same as for (CB4, CB5) 
to VB6 (but correspond to a higher radial excitation) , so 
they are both allowed in D 6h symmetry, and CB9 — VB7 
is a LH transition purely due to valence band mixing. 

E. Summary of optical transition spectrum 

Optical transitions with polarization perpendicular to 
the nanowire axis are clearly more dominant for the low- 
est energy states. This is expected due to the oblate 
(quantum- well-like) aspect ratio of the QD (Fig. 1), the 
other optical transitions seen with optical activity polar- 
ized along the nanowire axis are linked with valence band 
mixing and overlap between the conduction band and 
the LH components. Accordingly, we see an increasing 
number of transitions with this polarization for higher 
levels, as the valence band states become increasingly 
LH-like (see Tab. V). A computation with more energy 
levels, or a narrower QD with stronger lateral confine- 
ment, would have allowed to reach the ground LH-like 
level with strong dominance of z— polarized optical tran- 
sitions. Our approach would then evidence the second 
kind of DPGPS. If the QD aspect ratio would be re- 
versed to prolate a reversal of the roles of HH and LH is 
expected. 

Excitonic effects, which we have neglected in the 
present work, are not expected to change very much this 
global picture of polarization anisotropy. It would intro- 
duce electron-hole exchange effects which would split all 
the dominant transitions into a further observable fine 
structure (which can be seen as quadruplets due to spin 
degeneracy). The symmetries of the fine-structure exci- 
tons can be easily obtained by the product of irreps 31 . 
Sometimes doublets may remain, other times all lines 
would split in D 6 h- 

In a real situation, note that the dielectric mismatch 
between the high index nanowire and the surroundings 
may filter emission along the nanowire axis, altering the 
intrinsic optical response of the nanowire QD discussed 
here. 



VIII. CONCLUSION 

We have presented a systematic procedure based on 
class operators for symmetry analysis of the electronic 



states of a QD. The procedure, called PTCO, is based on 
postprocessing and alleviates the need for a code special- 
ized with respect to a given symmetry (which would how- 
ever feature significant gains in memory/time 37 ). PTCO 
has been demonstrated using the k-p method for the con- 
duction and valence band of hexagonal GaAs QD grown 
within Alo.3Gao.7As nanowires. The high D 6 h symme- 
try of the QD heterostructure is partly broken by the 
C3V crystal symmetry carried by the Luttinger Hamil- 
tonian, and by the unsymmetric computing grid. Using 
the PTCO on the computed results, we have been able 
to sort all these symmetry breaking effects, and quantify 
them. We have demonstrated that the numerical grid 
effect was small and could be compensated by PTCO, 
whilst the deviations from the approximate elevated D 6 h 
symmetry towards the true C3 V symmetry, albeit small, 
could be measurable. 

The PTCO is simple to program, intrinsically system- 
atic and automatized, and is carried out in a single step 
for a given symmetry, delivering at the same time sym- 
metrized states with corrected energies and a classifica- 
tion for all quantum states. It is flexible, applies inde- 
pendently of the method used to compute the electronic 
structure, and can be tuned with very little efforts to 
analyze a higher symmetry in a second run. This en- 
abled us to investigate the proximity of each state to an 
approximate elevated symmetry. 

In a second step, all the quantum states were an- 
alyzed using projection operators to give quantitative 
weights for every symmetry group, for every ultimately 
reduced envelope function present in every spinorial com- 
ponent. We could then verify all analytical predictions 
made concerning the UREFs in the elevated D 6 h sym- 
metry (Eq. (B.6)). This approach allowed to identify 
dominant DPGPS and dominant UREFs for all states 
considered. 

The analysis of the dominant UREFs opened the possi- 
bility to attribute additional azimuthal and radial quan- 
tum numbers to every state. The natural sequence al- 
lowed to explain the order of irreps in the computed re- 
sults, by subduction from 0(3) to D 6 h- 

A final summary of the classification is given for all 
valence band states in Tab. V, and for all conduction 
band states in Tab. VI. The classification provides in- 
sights into the origin of degeneracies and quasi degenera- 
cies, and allows to predict all selection rules and most of 
the polarization properties. The information can also be 
used to construct approximate product states, relevant 
for the interpretation of excitonic and multi-excitonic fine 
structure 31 . 

Finally we were able to interpret completely all the de- 
tails of the computed optical spectrum with the help of 
the classification of the states. We unveiled a large num- 
ber of missing transitions, which were shown to stem not 
only from approximate elevated symmetry but also from 
approximate azimuthal and radial quantum numbers. As 
valence band mixing was included in the model, we could 
identify all the corresponding fine structure, responsible 
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for some spread in the optical transition spectrum, and 
to verify analytical predictions concerning the appear- 
ance of "magic ratios" in polarization anisotropy. 

High symmetry GaAs nanowire QDs as investigated 
here are particularly interesting. First, since there are 
a large number of forbidden transitions in their optical 
spectrum, the multi-excitonic spectrum is also simplified 
and can be better interpreted. Second, the doubly de- 
generate exciton states due to symmetry are well suited 
for the generation of entangled photons. 
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Appendix A: Choice of standard irreducible matrix 
representations for C3 V and D6h 

For one dimensional irreps T, the set of matrices 
{D r (g),g £ Q} are simply the characters of the respec- 
tive irreps i.e. -D r (g), x r g, Vg, where x r (ff) is listed in 
26. When r is of higher dimension, the matrices {D r (g) 
are only unique up to a similarity transform U(T), i.e. 

D r (g) = U(T)-'[Dl lt (g)] H U(T) (A.l) 

where D^ lt (g) is a matrix irrep listed in 26, and U(T) is 

restricted to be unitary. For our purpose the [-D^itO?)] 
are not the most numerically convenient since they diag- 
onalize rotations instead of a mirror operation w.r.t. a 
symmetry plane, e.g. a v \ (Fig. 2). We thus choose the 
following to define our D T (g) matrices: 
C 3v , irrep E: 



i_ 1 

V2 V2 



(A.2) 



C3V, irrep E X j 2 (with case A in Ref. 26): 



V2 V? 

1 —i 
V2 V2 



(A.3) 



D 6h , irreps Ef. 

Uv 6ii (E i ) = U C3v (E), i = 1..2 (A.4) 
D 6h , irreps E jtk : 

1 5 

U^{Ej, k ) = U C3v {Ei/ 2 ), j = -,-,k = u,g (A.5) 



UD 6h (E 3/2tk ) - 



V2 sf% 

— 1 — i 

V2 V2 



, k = u,g 



(A.6) 



Symmetrized bases 26 , also called Symmetry Adapted 
Functions (SAF), follow from subduction from 0(3) to 
Dgh and allow to find easily the proper HSBF basis corre- 
sponding to the top (bottom) of the valence (conduction) 
band, respectively, of diamond or Zinc Blende. Unfortu- 
nately the SAF in Ref. 26 do not have regular properties 
under time-reversal, so some care is required. 

Let us first consider the bottom of the conduction 
band, and use linear combinations of zone center Bloch 
function denoted , m = ±|, corresponding to a 

"quantized axis" with z along [111] and x along [112] 
Fig. 2. They are nearly odd under spatial inversion in 
GaAs (hence the [...)* symbol as in Ref. 26). Subduc- 
tion from 0(3) tells us that the irrep E 1 / 2 u of D@h must 
be associated with this subspace. To construct the HSBF 
we can either use the SAF of Ref. 26 which are suitable 
in this case, and the change of basis corresponding to 
Eq. (A.5), or equivalently decide to diagonalize the set 
of Wigncr matrices corresponding to the symmetry op- 
erations (parametrization and factor system of Ref. 26), 
and find their reduction to the set D El / 2 - u {g) given in 
Fig. 19. The resulting conduction band HSBFs are: 




(B.l) 



For the top of the valence band we use the set of similar 
Bloch- functions labeled |§,ti) , m = — nearly even 
under inversion. Then the 0(3) subduction tables to Dgh 
indicate reduction to the E 3 / 2 g + E 1 / 2 g irreps. This time 
one cannot use all the SAFs of Ref. 26 if one wants to 
keep an invariant form of the time reversal operator in 
the HSBF, therefore we follow the route of reducing the 
Wigncr matrices corresponding to symmetry operations. 
Using a suitable, but freely chosen, set of phase factors 
compatible with the irreps of Fig. 19, one finds the fol- 
lowing D@h valence band HSBFs: 



\E- 



3/2,37 



1 \ _ i 13 3\ 1_ 13 _3\ 

V v/2 12'2/ v/2~l2> 2/ 



1^3/2,3, 2) - 75 



V2 I 2' 2 



3 3\ 1_ 

2' 2/ v /2 
3 3 \ _ i 
/ ,/2 



and 



1/2,3 
1/2,3 



.1) 

2) 




(B.2) 



(B.3) 
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This choice satisfies three desirable constraints: 1) the 
matrix form of the time reversal operator is invariant 
when the valence band basis is ordered as 

{\^3/2,g: l) ! \ E l/2,g, l) j |-^l/2,gi 2 ) , \E 3 / 2t g, 2)} (B.4) 

hence it preserves the canonical p, q, r, s form (Eq. (II. 3)) 
of the Luttinger Hamiltonian in this basis, 2) it is also si- 
multaneously a HSBF for C 3v , 3) their transformation 
laws are given by our set of standard representations 
listed in Fig. 19. 

For C3V the choice Eqs. (B.2) and (B.3) is in agreement 
with Ref. 24 with the correspondence 



\E 3/2ig ,i)^\ i E 3/2 ),i = 1..2 
whilst it differs from Ref. 1 by 



(B.5) 



E 3/2,g, i 

Ei/ 2t9 ,i 



(-iy 

(-!)< 



3-i 



E3/2) 



'3/2, 



,i = l. .2. (B.6) 



Although the matrix representation for C3V in Fig. 19 
are the same as in Ref. 1, we had chosen an opposite 
projective factor system for the improper operations in 
Ref. 1. The present projective factor system is now the 
same as in Ref. 26. 



Subduction rules to C3V (according to the case A in 
Ref. 26) will clearly reduce the spinors Eq. (C.l) to 
Eq. (IV. 3). Since Ei, E 2 — > E one might wonder why 
the partner function indices are reversed for E\ w.r.t. 
E 2 (or E in C 3v ) in Eq. (C.l) (Eq. (IV.3) in the main 
text). The reason is simple and lies in our choice of ma- 
trix representations in Fig. 19, where the representative 
of <j v i has a different sign for E\, reversing the a v i parity 
characteristics. 

Finally time reversal symmetry will bring further con- 



straints on the UREFs, e.g. via K^_ x 



1/2,11 



and 



-^j 1/2,11 due to Kramers degeneracy. Prop- 
erties similar to Eq. (IV. 4) can then be readily obtained, 
but are omitted for clarity. 



V2 



/ <t>f^ 

02 



1 

71 



4*2 



02 
*2 " 

v 

A E ^ 
01 



l_ A E ^ 
+ $f V 



\ 



- $ 



b 2 



(C.la) 



(C.lb) 



Appendix C: UREF decompositions of the valence 
band spinors in the elevated D6h symmetry 



ti 



3/2, it _ 







A 1." - B 2." \ 


1 




$f 1 -» - $f 2 " 


71 













(C.lc) 



Using the HSBF derived in Eq. (A.6), Eqs. (B.2) 
to (B.4), it is possible to decompose each spinor into 
UREFs with the help of Eq. (48) of Ref. 1 and the 
Clebsch- Gordon coefficients linked with the chosen stan- 
dard matrix representation of Fig. 19, and involving a 
minimum of arbitrary phase factors. In case of C3V 
symmetry one finds that the UREF decomposition of 
the main text Eqs. (IV.3), which remain identical to 
Refs. 1,24. In case of Dgh symmetry, and in the HSBF ba- 
sis order given by Eq. (B.4) below, one obtains the unger- 
ade valence band spinors given by Eq. (C.l). The corre- 
sponding gerade spinors have similar expressions where 
u — >• g. To lighten the notation we have omitted the main 
spinor index and the HSBF index in the UREFs, but they 
can be retrieved easily from the main spinor and the po- 
sition of the UREF (using Eq. (B.4)). For clarity the 0/<I> 
UREFs are not the same functions in Eqs. (C.lc)-(C.ld) 
w.r.t. Eqs. (C.la)-(C.lb), and we have used the capital $ 
to distinguish the E x and E 2 UREFs in the E 1/2>g HSBF 
components from the £3/2,3 HSBF components. 
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71 



x pE 5/2u 
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71 
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A E 1, 



A E L 
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